SOME APPLICATIONS OF FIXED POINT 
THEOREMS IN CONTROL THEORY 


A Thesis Submitted 

In Partial Fulfilment of the Requirements 
for the Degree of 

DOCTOR OF PHILOSOPHY 


BY 

SUBHENDU DAS 


- » to the 

DEPARTMENT OF ELECTRICAL ENGINEERING 

INDIAN INSTITUTE OF TECHNOLOGY KANPUR 

OCTOBER 1972 



OERTIFIOiiirE 


Certified that this 7/orlc , ’ • Some Applications 
of Fixed Point Theorems in Control Theory'hy 

Mr» Subhendu Das, has been carried out under . 

% 

my supervision and that this v/ork has not' - 
been submitted elsewhere for a degree . 



(Dr. B. Sarkar) 
iissistant Professor 
Department of Electrical Engineering, 
Indian Institute of Technology, 

, • Kanpur. 


I SI office 

I for th approved 

I Decree of 

I Doctor of Philosophy (PJi.D.) 

I accordance with the 

f regulations of the Indian 
^mmute of Technology ,„a„p„ 
Dated; j .15 





3.b' 


\5^‘l- D- \)PS-SDtv\ 



ACOOWLEDCTIEM S 


The author is grateful to Dr. B. Sarkar, 
Assistant Professor, Department of Ele- 
ctrical Engineering for allowing to work 
in an atmosphere of dignity and confide- 
nce and extending various kinds of help 
during his stay here. He is also profo- 
undly indebted to Dr. P.O. Das, Assistant 
Professor, Department of Mathematics for 
many valuable discussions and suggestions 
on this subject which led him outside the 
traditional realm and view point of engi- 
neering. Thanks are also due to Mr, S.S, 
Pethkar for typing the thesis. 



SYNOPSIS 

SUBHENDU Ph.D. 

Indian Institute of leclinology 
M'PUE 

October, 1972 

SO&IE ^IPPLICIIDIOIS OE EHEB POINT THEOREMS IN 
CONTROL THEORY . 

Many control problems can Be represented in tiie form 
of Integral Equations (IB). Noting tliat these are special 
cases of operator equations, a functional analytic view 
point (mainly the constructive approaches) to study the 
qualitative properties, like existence and uniqueness of 
solutions, of these equations, has been adopted. Emphasis 
has been given to show hov; these techniques can be applied 
to control problems ra,ther than giving best possible 
results or solving complicated large practice oriented 
problems. Several numerical . examples have been included 
to demonstrate the procedures. 

For Boundary Yalue Problems (BVP), of both t?ro point 
and multipoint cases, a simple equivalent IE represen- 
tation has been given. The technique involved can 
accomodate large class of constraints on the system 
dependent variable (called state) and other parameters. 
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Usirig this technique, as illustration, a yery general nour-' 
linear controllability problera has been given an equi- 
valent IE representation, ffith the help of some modifi- 
cations of Banach* s fixed point theorem, the equivalent 
IE for the BVP has been studied. 

The most important aspect of this work is the study 
of the controllability problems. Eor this problem, our 
approach is motivated by the fact that these are essentia- 
lly two point BVPs, where an unkno^m parameter, called 
control, is to be found out, ?;hich will steer the system 
between any given pair of states. A set of lEs has been 
given, the solution of which represents a solution of the 
controllability problem. Using fixed point theorems on 
these IBs existence of atleast one control has been assumed 
The metliod presented here also finds out the numerical 
values of the control, A very simple result has been 
obtained for the controllability of non-linear system, 
when the dimensions of both the state and control vectors 
are identical. 

The stability problems have also been investigated by 
fixed point theorems. Eor the non linear forced system, 
when the forcing function independent on both state and 
independent variables (called time) of the system, Banach’s 



fixed point theorem has been nsed. hnd for the open 
loop systems, that is, v/hen the forcing function is 
dependent only on time, a modification of Kalcutani’s 
theorem to Banach spaces, given by Bhonenblust and 
Karlin, has been used vathout any convexity conditions 


on the system 
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In this chiapter we present the necessary 
mathematical hackgrotoid required for the understanding 
of the following chapters. The reader wishing to obtain 
a more detailed of these concepts should refer to [^1 - s] 
We will be using the following basic set theoretic 
notations : ■ • 


implies 

implies and implied by 
G belongs to 
Cl subset of 


^ such that 
a there exists 
y for all 

iff if .' and only if 


BMAOH SPACES 1.1 
Linear Spaces 

Let X be a non-empty set. Assume that 
(i) there is a rule that v/ill allows uS to construct, 
for every two elements x,y S X another element 
2 e X, called the sum of the elements x and y and 


denoted by z = x + y- 
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(ii) there is a rule that will allovr us to construct for 
every element x G X and every scalar a an element 
u G X called the product of the element x and the 
scalar a denoted hy an = u. 

DEFIhIT IO I'I ~ l.l.l 

A set X is called a linear space if 
additions and scalar multiplications are defined 
; on X and the follov/ing rules hold; 

(1) x+y = y + x for all x, y G X 

(2) (x+y)+z = x+(y+i 2 ) for all x,y,z G X 

(3) There exists an element 0 G X (the aero element) 
such that X + 0 = X for any x G X. 

(4) Eor any xGX, SyGX 9x + y = 0. 

(5) lx = X for all x G X. 

(6) a(px) = (ap)x for any x G X and any a and p. 

(7) (a+p)x = ax + px for any x G X and any a and p. 

(8) a(x+y) = ax + ay for any x,y G X and any a. 

A linear space;' is called a real or 
complex according to whether the scalars are the real or 


complex number system 
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DEFIIIT ION 1.1.2 

A linear space is called a nomied linear 
space if a rule exists, which, associates v/ith 
every element x G X a real number (called the norm 
of an element x and denoted hy [|x||). 1 his rule 

must obey the' folio v/ing ; conditions (norm axioms) ; 

( 1) 11-11 >0 , 1 1 x| I =0 '4=^ X = 0 . 

(2) llx+yll < Mxll + llyll- 

(3) 11 axil = 1 aldlxlD 

where 9C is any scalar. 

Henceforth we will asstime that Z is a 
normed linear space. 

Open _and_ Closed. Sets: 

hSFIN IllON 1.1.3 

If Xp is a point of Z and r a positive 
real number, then the open sphere S^(Xq) with 
centre x^ and radius r is the subset of Z 
defined by 

S^(Xq) = I X : llx - XqI 1 < r}; 
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dekhttioe; 1.1.4 

jj. sul'>se't G of X is calle-d an open set 
if given anjr point x in G there exists a 
positive real number r such that S 2 ,(x) c;;;; G i.e. 
if each point of G is the centre of some •pen 
sphere contained in G. * 


TliEO RM 1.1.5 

In any normed space, each open sphere is 
an open set . 


TJ-IEORffi 1.1.6 

A subset G of X is open it is a 
union of open spheres. 

TifflOM 1.1.7 

jiiiy \inion of open sets in X is open and 
any finite intersection of open sets in X is open. 

DERIHITipi 1.1.8 

Let A be a subset of X a point x in X 
is called a limit point of A if each open sphere 
centered on x contains atlea,st one point of A 


different from x 
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LEI'IMTXOH - 1.1,9 

E subset P of X is called a closed set 
if it contains eacii of its limit points . ' 


A subset P of X 
ment p'^ is open. 


Tpogai 1.1.10 

s closed its compie- 

THEORai 1.1.11 ■ 


Any intersection of closed sets and any 
finite union of closed sets in X are closed sets. 

DEPnn T 10 1 1.1.12 

If is a point in X, and r a non- 

negative real number , then the closed sphere 

with centre x^ and radius r is the subset of 
o 

X, defined by 


Sj,[xo] =' {X : 1 [x - x^[ [ < rp 


THEORIM 1.1.15 


In a normed space each closed sphere is 
a closed set. 


Gompl eteness : 


Iiet Xj^ = { x^ f ^2* * * ‘ j • .j 

be a sequence of points in X. We say that {x^> is 
convergent if there exists a point x in X and for 
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each G > 0 there exists a positive integer such 
that 

^ > ^0 I 1^11“ ^li i e . 

Every convergent sequence {30^^} has the 
following property: for each G > 0, there exists a 
positive integer n^ such that in,n > n^^^ | H 6. 

A sequence with this property is called a Cauchy sequence . 
It is to be noted here that not every Cauchy sequence is 
convergent because for convergence it is necessary that 
the limiting element should belong to the set. 

hEEIIlI TIOH 1.1.14 
A normed space is complete, if every 
Cauchy sequence in it is convergent. 

TIiSpRSI 1.1.15 

Let X be a complete normed space and let 
Y be a subspace of X. Then Y is complete it 

is -closed. 

pEraiTIpI ■' 1.1.16 
A Banach space is a complete normed 
linear space . 

;ife ?H1 use the ?rord ’B-space’ and 
•Banach space' int e r change bly. 

pPERATORS: 1.2 

In what follo?ra, ?/e assume that X and Y 


are Banach space s . 
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pEgl llTIOH 1.2.1 

- 11 . function A is callod an operator if to 
eacli elenent x in some sutset D of A assigns 
one element y G I . The set D is called the domain 
of A, while the- set of all elements y = fe for all 
X G D is called the range of A. 

DEfimiPI 1.2.2 

An operator A is called hounded in X if 
there exists a constant M such that | | Ax| [ < M| |x[ j 
for all X G X. 

fEPIl'IITIpI 1.2.3 

Am operator A is called continuous if for 
any G > 0 there exists a 6 > 0 such that jjx^-Xgij 

< 6 implies ^ ^ 

pEEIIIITIpl 1.2.4 

An operator A is called additive if 
i-Cxj^+Xg) = Aoi^+mcg every x^^Xg G X, and 
homogeneous if A(a3c) = aAzz for everj^ scalar a 
and every x G X. jln operator v/hich is hoth 
additive and homogeneous is called linesn, 

DEEIhlTipi ■ 1.2.5 

The norm of a bounded linear operator A 
is the smallest number M which satisfies [|Aoc[ [ 

< M 11x11 for all X G X. It is denoted by 1 1 Al I - 
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It can be sliown tbat j jA] | is given by 

1 jAxl 1 

p === Sup II Axil 

1 11x11=1 

where Sup stands for supremuni i^e. the lea,st upper bound. 
Following theorems give some important properties of 
linear operators in Banach Spaces. 

lIlBOIffiM 1.2.6' 

If A is a linear operator from X into Y 
which is continuous at some x^ G 1, then A is continuous 
at every point x G X. 

I HE ORM 1,2.7 

Let XjY be Banach spaces and A a, linear 
operator from X into T. Then A is bounded ^^A is 
continuous . 

T HEO REM 1.2,8 

Additive and continuous operators are 

linear. 

^ [X>Y] we will denote the class of all linear conti- 
nuous operators mapping X into Y. [_2,Xj will be denoted 
by [X] . 

^pRIM •; . 1.2.9 

Let X and Y be Banach Spaces. Let the 
norm of A G [X^Y] be defined as in Definition 1.2.5. 




Then the space [Z,Yj Is a Banach space with hespect to 
the point wise linear operations. 

DEBlhlTIOII 1.2.10 

let A he an operator mapping Z into Y, 

¥e say that the operator is an onto mapping if corres- 
ponding to any element y 6 Y E-; an x 8 Z 5 Ax = y. If 
under the operation A two different elements in Z 
always have different images then A is called an one- 
to-one mapping of Z into Y. 

TH EOREM 1.2. 11 

Let Z and Y he Banach spaces and A a 
linear operator mapping Z onto Y. Let a positive 
nuuaher-mii- .exist such that for every x e Y 

1 1^1 1 > 1 

Then the operator A has a linear inverse A , where 

MEASTOE AM D IMBSEmiOH l.S 
Following is a brief exposition. of the 
modem theory of integration. We follow the lines of 
Halmos [5] . 

Measurable sets 

DEFIHITIOH 1.3.1 , 
Let Z be any non empty set. A family S of 



10 


subsets of X is said to be a a -a.lgebra if it satisfies 
the followiiig conditions 

(i) 0(empty set) and X e S 

(ii) Aes, BeS=^. A-BeS 

(iii) if fA^> is an infinite sequence of sets in S then 

U 

n = 1 


We assinne S to be the smallest cy - algebra on X, 

BBBllIlIOh 1.3,2 

A measure is an extended real valued set 


function ^ , defined on a o -algebra S such that 


(i) V (E) > 0 E e S 

(ii) y (0) = 0 

(iii) S, A^f) \ = 0 

for m n. 


oo 


oo 


li (U^)^ 2>(e ) 


n=:l 


n=l 


E EEiroiOB 1.3.3 

An outer measure u * for any arbitrary set 
E in S is defined by 

v*(E) = Inf {mCP): p=>b, E G S } 

In a similar way an iimer measure Hjj. is defined by 

d^(E) = Sup {f'(E): ICE, S’ 6 S } 

DEEINITIOI 1.3.4: 

A set E in S is said to. be measurable if 
its inner and outer measures coincide with its measure 
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i.e. if y*(B) = y (E) = y ^(S) 

It call Toe siio’wn. that all the sets in S is measurable, 

■ . A statement is said to hold almost everywhere (a.e) 
if it is" true:: iexGBpt over a set of measure zero. 

PBglllll Ol 1.3.5 

A sequence {•’ f^} of a.e. finite valued 
measurable functions converges in measure to the mea- 
surable function f if for every e > O" . .. 

lim w( X : If^(x) - f (x) [ > e ) = d. 

n ^ 

lEPl hlllOIif 1 .3 . 6 

A measurable space is a, set X equipped with 
a a - algebra S of subsets of Z . 

let Z be the real line, P, the class . 
of all bounded semiclosed interval of the form [a,b ), S 
the cr _ algebra generated by P. let be the set 
function defined on the sets of S by Vta ,b) ) = b - a. 

Then it can be shown that y is a measure on S, This 
measure is known as lebesgue measure. Obviously this is 
a particular case of measure, so most of the results hold 
for this case, without probably any change. 

I ntegra ble Punction 

DBPINITION 1.3.? 

A real valued function is measurable if 
for any real number c the set {x: f(x)<c } is measurable. 
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DSmi TIOU 1.3.8 

A function f defined on a measurable space 
X, is called simple if tbere is a finite disjoint class 
{E^, . . .E^}’Gf measurable sets and a finite set 
{ • " -a^} of real numbers such that 

f if m G i = l,2,..,n 

= I 0 if xf; Bg.... 

I HEOREla 1.3.9 

Every extended real valued measurable 
function f is the limit of a sequence { f^ } of simple 
functions . 

E EEIEIIION 1.3.10 

A simple function is integrable if u (Bj^)< °° 
for every index i for which 7 ^ o. The integral of 
f is defined by 

H 

/ f d Ti= 2 a. 11 (E.) 
i=l 

DBEII mITI OI 1.3.11 

A sequence of integrable simple functions 
is said to be mean fundamental if 

y (f^, fj '-* o as n, m - 
where u (f jg) = /| f-gj b.y 
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DEraiTIOil 1*3*12 
iiii almos-fc everywhere finite valued 
measurable function is said to be integra^ble if there 
exists a mean fundamental sequence { f^^} of integrable 
simple functions which converges in measure to f * 

The integral of f is defined by 

/ f d6:,. = limi / f du , ' 

H -*■ oo ^ 

TIEORM 1.3.13 

A bounded measurable function is integrable, 

, DIEFBRE lT Ii l E''yj.aiQirS 1.4 

In this section we present some standard 
results on initial value problem of ordinary differential 
equations which v/ill be helpful in understanding boundary 
value problems. 

L ipso h itz Gondition 
Let f be a function defined over the 
set SC R ^5 v/here R^ denotes the real n dimensional 
space,v/ith values in R^. Y/e say that f satisfies the 
Lipschitz (Lip) condition on S if there exists a constant 
L > o such that 

l|£(x) - f(y)ll <111^ - ytl 

for all i'x,y e S . The constant L is known as Lipschitz 
constant. A function satisfying Lipschitz condition is 
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called lipscliitsian. I'lote that a Lipschitaian function 
is continuous . 

If the function f posseses first partial 
derivative v/ith respect to all the components of x over 
its domain of definition, then it can he shown, that the 
lipscliitz constant L can he tsilcen as 


n 

I = max E Sup | a. . (x)j 
i 0=1 X 


?7here a. . s "are components of the Jacohian matrix of f. 


Standa rd Res ults 
We assume the following properties for 


the function f . 

(i) f: I X I) - R ?7here I and D are open subsets of R 
and R^ respectively. 

(ii) f is continuous in x G D, for each fixed t G I and 
meanurahle in t G I for each fixed x G D . 

(iii) there exists a Lehesgue integrahle function ‘ m(t) 
such that i[f(t,x)j[ < m(t) for all x G D and almost 
everyv/here in I. 

The conditions (ii) and (iii) are known 
as Carat heodory conditions. 


DEPIIIT IO R 1.4.1 
let S he some subset of I x D . Then by 


the solution of the initial value problem 
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± = f(t,x) x(to) = a (1.4.2) 

v/e mean an alDsolutely continuous function 0 ( 1 ) such 
that 

(i) (t,0(t)) e s , 

(ii) = a 

(iii) j^(t) = f(t, 0 ( 5 )) foi a-11 t , except over a 
set of Lehesgue measure zero, of the interva,! for 
v/hich 0 ( 5 ) is defined. 

LmU 1.4.5" 

t 

The initial value problem (1,4.2) is 
equivalent to the following integral equation 

t 

x(t) = a + / f(a,x(s))ds (1.4,4) 

The above lemma a,ssures that instead of 

studying the system (1.4.2) we can study the integral 
equation (1.4.4) which in some cases may be easier to 
handle with. 

TlffiORM 1.4.5 

let f(t,x) satisfy all the conditions 
mentioned above. Then there exists a solution 0 ( 5 ), 
defined over some interva,! of I, of the initial 
value problem (1.4.2). In addition if f(t,x) 
satisfies the lipschltz condition in x for all t e l, 
the solution is unique. 
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JilLUB PRp BLMl 1 . 5 

Tile linear boundary value prolalem for 
ordinary differential equation has been studied almost 
exhaustively. Any standard hook on differential equa- 
tion will give good record of these results, see for 
example [6,7j . The material presented in this section 
closely follows falh [8] . 

Consider the linear two point hound ory 
value problem (TPBVP) • 

x(t) = A(t) x(t) + f(t) (1.5.1) 

Mx(o)+ Nx(l) = a (1.5.2) 

where x(t), f(t), a are n-vectors, A, M and H are nxn 
matrices. Assume that A(t) and f(t) are defined and 
mea,surahle over an open interval I containing the closed 
interval [o,l')], also that there exists an integrahle 
function m(t) on I such that | |A(t)l j < m(t), 

1 |f(-fc)! 1 < m(t). 

DEP IlIT IOh 1.5.3 

An n- vector valued function 0(t) is 
called a solution of the above TPBVP if 

(i) 0(t) is an absolutely continuous function 

(ii) 0(t) == A(t) 0(t) + f(t) for almost all t in[0,l] 

(iii) M 0(o) + B0(1) = a . 
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let ^(t,t^) le the Hmdamental matrix of 
the linear system k = A(t) x with $(tQ,tQ) = I. let 
us assume for coiivinience that to = o and t G 
Then vm know that 

t 

x(t) = «>(t,o) x(o)+ / < 5 (t }.T)f( -r)dT (1.5.4) 

0 

is the solution of the differential eqr.ation (1.5.1) for 
some given initial condition x(o). If this solution isiso 
satisfy the Bounda,ry condition (1.5.2) then we must have 

1 

Mx(o) + h[ f( 1, 0 )x(o)+ / f(l,T) f (t )dT ] = a 

o 

1 

or [M+H $(1,0)] x(o) = a- N / ^(lj,'T)f ( T)dT = b(say) .. 

o 

Thus Yfe have the follovfing lemmas . 

ISm 1.5.5 

let S denote the set of all solutions of 
the TPBVP mentioned above. Then S is non empty iff 

1 

a - K / ^1, t) f( T)dT =5 
0 

is an element of the range of the operator (M+I <e(1,o)] . 

lEffilA 1,5.6 

If the matrix [m +11 $(l,o)] is non singular 
then the IPBA/P defined above has an Yinique solution 
viThich can be expressed in the following form. 
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1 

x(t) = H(t) a + / G(t, s) f(s) ds (1.5.7) 

o 


v/here the Green’s Matrices H(t) and G(t,s) are given "by 

-1 

H(t) = ^ (t,o) [M -t- h >$» (l,o)] (1.5.8) 

( H(t) M ^(o,s) 0 ■< s < t 

&(^^»s) =< (1.5.9) 

I -H(t) hf (l,s) < s < 1 


respectively. 
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CHiPlM - TWO 

FIXED POUT THBORS![S 

I RTRODUGTIOF 2.1 

jto. operator A mapping a Banach space 2 

into itself is said to have a fized point if there exists ■ 

^ * 

an element z G Z such that 

Az = z 

i.e, the element z* remains invariant under the operation 
A. The study of the fized point of an operator is eq_ui- 
valent to the study of the solution of the operator 
equation 

Ax = z (2,1.1) 

or equivalently 

Bx = 0 (2.1.2) 

where B = A - I and 0 is the null element of the B-space 

Z. 

There are various theorems, knovm as 
fixed point theorems, available for the study of the 
existence and uniqueness properties of the fixed points 
of an operator. Most commonly used results of them are 
the Schauder’s [9] and Banaches [9] fized point theorems. 
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ITevrtons metliod 1 ] or the Bellman’s quasilinearization 
technique [ 10 ] has also "been widely used [ 11 J ior 
solving the operator equation (2,1,2). It has heen shovai 
that this method is a special case of Kantorovitch’ s fixed 
point theorem [ 1 ] . 

In Schauder’s fixed point theorem one 
needs to find a convex closed compact set which will he 
mapped into itself hy the operator. This free choice of 
convex set mal-ces the theorem very ■ flexible. But in most 
cases it becomes almost impossible to find out a siiitable 
convex set, specially if the operator is not of general 
standard form. Moreover this theorem does not give, any 
constructive procedure to compute the fixed point, which 
is quite often required in applied analysis. 

On the other hand Banach’ s theorem, 
popularlj’’ laiovm as Contraction Mapping Principle (CMP) 
is very simple to apply and it gives a constructive 
procedure to find out the fixed point, Ihe method , 
generates a sequence of elements which converges to the 
soliition. The theorem requires that the operator should 
be Lipschitzian [ see section l.i ] with lipschitz 
constant less than one. It has been found that many 
operators exhibit solutions even though their Lipschitz 
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cons'fcan'ts are not less than one. This means that the 
Banach’s theoiem usually puts a Yerj strong condition on 
the operator. 

To avoid this difficulty many modifications 
have heen presented. In this chapter vre present some of 
these results, virhich are helpful in solving problems. A 
very interesting result in this connection is the one due 
to L.M. Grave [ 12 ] . We shoY/ here that the modified 

contraction mapping principle as given "bj Palb [ 8 j is 
nothing but a special case of the Grave's theorem. Y/e 
will be using these results for the analysis of Boundary 
value, Qontro liability and Stability problem.s for the non 
linear ordinary differential equation in the following 
chapters. 

The other widely used method, the Bevrbon’s 
method, has a fa.ster rate of convergence of the sequence 
it generates, than that of CJIP, but at the expense of 
large number of computations, specially in the evaluation 
of the deriva.tives, an oftenly unwanted computation. 
Otherwise tliis method does not have any special advantage 
over the Clip. Indeed the problems that can be solved by 
CldP can not alv/ays be solved by Newtons method, on the 
other hand the converse is true, A detailed and comparative 
study of these methods can be found in [ 8 ] . 



BANACH'S THEOREM 2.2 


Throughout in this chapter v/e assime that 
X is a Banach space (B-space) with the norm represented 
hy jj .j| . By small letters x,y, arbitrary general 
elements of X will be represented. Capital letters A 
and B will denote .-eperators defined ever X with values 
in X. 

THEOREM 2.2.1 

(Banach) 

let A:X -* X be an operator satisfying the conditions 

(i) ASCS where S is a closed subset of X 

(ii) j I Ax Ayj \ < a 1 jx - yj | for all x,y e S and 

with o < a < 1 


Then the operator A has a unique fixed point in S i.e, 
there exists a unique element x G S such that 

Ax = X* 

And based on any x ■ G S the sequence generated by 


“ ^n 

converges to x* with respect to the normed topology 
of X. 

In many cases, because of the method used, 
it turns out that the condition (i) of the above theorem 
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puts a stronger restriction on a than virhat is mentioned in 
condition (ii) specially v/hen the operator is not defined 
over the whole space. 

Suppose S = Sj, a closed sphere of radius 
r with centre on the origin, and let A he defined over 
with AD = a. Then for any x 8 S condition (i) requires 
that Ax G S . We have 

I 1 Ax| { = 1 I Ax - AO + aj I < j j iix - AO j | + j j a| ] 

using condition (ii) in the first term of the right hand 
side of the above esspression, we get 

1 I Axl 1 < a| jx - Oj 1 + 1 1 a| I < + | ] aj | 

Thus for liX to belong to we require 

ar + 1 1 a| 1 < r 

Hall 

i.e. 

The above shows that a should be usually less than what 
is demanded by condition (ii) . lu vie v; of the above 
fact in most of the cases it is sufficient to verify the 
condition (i) . 

In future we will refer the above theorem 
as Banach’ s fixed point theorem or the contraction mapping 
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principle . The condition (i) is hnown. as self mapping 
property and the (ii) as contraction property. The 
condition (ii) is also kno?m as lips chit z condition with 
a as lipschitz constant. 

MODIPIOiTIOlS 2,3 

Is mentioned "before the restriction on a 
in GI\!EP is quite strong, following theorem due to Chu 
and Liaz [ 13 | and Kantorovitch [ 1 j gives a modifica- 
tion over the above result. 

TH EOREM 2.3.1 

(Ghu, liaz, Kantorovitch) 

Let A and B "be two operators defined over a non-. 

empty set S with values in S. Let B possess a, 

—1 

right inverse B . Then the operator A has a fixed 
point iff B ilB has a fixed point. 

The essential idea "behind the a"bove 

said thee rem' is that if A is not a contraction operator 

—1 

and if ■.'-there exists a B operator such that B AB is 

contraction, then the fixed point of A can "be computed 

-1 

using the CMP over the composite operator B -AB. It is 
to "be noted that if x* is a fixed point of B-^AB then 
Bx* is the fixed point of A. The. above follows from 

X* = B-^ilEx'" ABx* = Bx* 
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In some cases it majj' be possible to reduce 
tbe lipscliitz constant of tiie operator A. by suitably 
redefining the norm of the underlying space, so that the 
operator A becomes a contraction Y/ith respect to the new 
norm. It has been shown in [ 13 ] that this method is 
nothing but a special case of the theorem 2.3.1. 

Following theorem can also be considered 
as a modification over GWB. A very interesting applica - 
tion of this result can be found in 14 ^ * 

IHE OEEM 2.3.2 

(M)' 

The fixed point of the operator A is the fixed 
point of the operator and vice versa, for any 
positive integer n > 1. 

Thus if for some integer n the operator 
is contraction, then using CUP on , the fixed point 
of A can be numerically computed, 

GR AYED S THEOSa i 2.4 

The theorem due to L.M. Grave [ 12 ] may 
be considered as a significant modification over the 
Banach's theorem. Many problems which cannot be solved 
directly by using Banach's theorem can be solved by 
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G-rave*s theorem. This theorem is "based on the following 
lemma due to Banach [ 15 ] . 

imaiA 2.4.1, 

( ' Banach ) 

If B is a continuous linear operator of the space X 
onto the v/hole space Y, there exists a num"ber p > o 
such th8,t for each y G Y there is an x G X with 
y = Bx and 

11^11 < ^ llyll 

T HEORM 2.4.2 

(Graves) 

Let A "be a continuous opersvtor defined for j |3c| r 
ViTith values in Y and iiO = 0. Let B be another 
linean continuous transformation of X onto the whole 
space Y such that 

j j Axg - B(x^-Xg)} 1 1 a I ~ Xgj | 

whenever x^ and Xg belong:.- to S^, where l/a is greater 
than the constant ^ of lemma 2,4,1, Then the equation 
y = Ax has a solution x in whenever 

1 1y1 1 1 • 

It has been shown [ 12 ] in the proof that 
the sequence { } generated by the alternate application 
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of the formulas 

^n-1 ~ (2.4,3) 

- ^ + “^1 (8.4.4) 

where Xq = 0 and = y, converges to the solution of 
the equation 

y = jIx . 

The difficulty with the above theorem 
is to find out the value of , v/hen aiid ^ 

are known, by using equation (2,4.3), But if B is 
invertible this will be greatly simplified. In this 
case the value for p, of Banach's lemma can be talcen as 

p = II B-l|l 

If B is invertible the G-rave's theorem can be rewritten 
in the following form. Bor convenience let us use the 
notation 

A e (lip, S,a)<::^l \ d*x ^iyj 1 < a 1 jx - yj | 
for all X, y G S. 

THEORM 2.4.5 

let Z be a Banach space and S^, a closed subset 
of it. Let A i - Z end let B be a linear 
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invertible operator: defined over X with, values 
1 

in X. Let | jB"* | P* Assume that the following 
conditions hold 

(i) [a 4- B] G ( lip, S^j a) 

(ii) [a + B] Sj^d 

(iii) a§ < 1* 

Then there exists a unique element x G such that 

* 

4*x = 0 

and based on any Xq G-S^^ ;the sequence' gexasratPh by 

Xn = ^ (2.4.6) 

^n = S 

converges to x . 

Proof (Sketch); v 

Let Xq, x^,... belong to S^. Then 

= llB-lnll < ^IPnll = x^l 1 

< p. | - r 0 

1 I = 1 IB-Vn - B-bn.il 1 < Pi lyn- yn-ll I 

= Pl I [•^B] Xjj - [A + B]Xjj^j_ 1 1 < ap| l^'•^u_ll I 


Since < 1, 
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the above implies jis Gauchy. let ^ x*, 
obviously x Q because is a closed set. 

Define y =: [A +B] x . Because of the continuity of 
[iU-B] , which follows from condition (i). 


lim y^ — lim ^A+bJx^^ — j^Ari-B^ lim Xj^= [^AtB'jx =y 
and hence from (2.4.6) and (2.4.7) we have 

* r 1 * 

y = I A+Bl X i „ 

^ / . f ■=!;> A X = 0 

y = B X j 


Q.E.D . 

It is easy to see that the sequence| x^ | 
generated by (2.4.6) and (2.4.7) can be generated in a 
simpler way by 

^ ® ] ^n (2.4.8) 


Thus we have the follovnng theorem which 
Ba-lb calls as modified contraction mapping principle [8] 
and is also available in [ 16 } . 

T HE OREM 2.4.9 

(Ea,lb, Krasanoselskii) 

let A and B be maps of 2 into X. let B be 
invertible . Then the fixed point of B“1[A + B] 
is the solution of 


ilx = 0 . 
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Thus if B is linear, the theorem 2.4.9 
can be considered as a special of Grave’s theorem. In 
the follov/ing chapters we will use Graves theorem in the 
form given by theorem 2,4,5 or 2.4,9, depending on the 
convenience. We observe t heat these two theorems are 
identical when B is linear, 

E3La£P LE 2.5 

We illustrate the applica.tion of the 
above results by on example from integral equa,tions. 
Consider the space 0, the class of all continuous func- 
tions defined over the real interval [o,l] with values 
in Z, an arbitrary Banach space with norm j i •! 
the nom in 0 be defined by 

jlxjl^ = sup { 1 |x(t)| 1^: x(t) G 1, X e 0 } 

te[o,i] 

It can be shown that with the above norm the class 0 is 
a Banach spa.ce. 

Let the function B(t,x); [o,l] x S^— X, 
S^ciS satisfy the following conditions 

(i) fCtjx) is continuous in x G for each fixed , 
t G [0,l] and measurable in t G [0,l] for each fixed 
X G . 



(ii) Tliero exisi:s an in"begrable fu.iici:ioii 111 ( 1 }) such, 
that I lP(t,x)i 1^ < m(t) for all X e . 

^(tjx) Soitisfies th.e follo?/iiig uniform lipschit 
condition 

1 |5'(t,x)- P(t,y)| 1^ < Ij jx-y] for all x,y e 
and oai t G [0,l] . 

(iv) P(t,o) = 0 for all t G [0,l] . 

lot q(t) be a continuous function belonging to th.e , 
class C v/ith the norm 

11 1 1 Ic 1 ^ 

where a is any non-negative real number, ■ 

We a.re interested in studying the 
existence and uniqueness properties of the solutions 
of the following integral equation 

t 

x(t) = q(t) - x(l) + /P (s,x(s)) ds t e [0,l] 

o 

By solution we mean a function 0 G O 
viThich satisfies the above equation. In above x(l) is 
not knovai. Clearly because of the presence of x(l) 
direct application of Om will fail, however small 
the lip schitz constant (for P) may be . Pormally we 
have the follo?/ing result. 
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Assume that the function I'(t,x) and q.(t) sa.tisfy 
the conditions mentioned above. Let o < L < § - 
Then the above equa,tion ha.s a unique solution. 

To verify this we proceed in the 
following way ^ along the lines of t he theorem 2,4,5. 
Define the operator A by 

t 

y = Ax y(t) = q(t) - x(t) - x(l) -f / T’(s,x(s))ds 

0 

The obvious choice for B is 
y = Bx y(t) = x(t) + x(l) 

and hence [A + B] takes the form 

t 

y = [a + B] X y(t) = q(t) + / B(s,x(s))ds 

o’ 

Clearly B is linear and invertible, B’ 

is given by 

X = B”V x(t) = y(t) - I y (l) 

which gives 

Ihllc i 1 iblU = I liyiio 

thus 

p, = 3/2 , 

we have 


Dor any x,y G 0, y(t.)x(t) G 


f 
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1 1 [A + B]z (t) - [^l+B] y(t)l I . 
t 

= II /:-t'5'(s,x(s)) ~ B(s,y(s))]ds I 
t 

< / j|B(s,z(s)) ~ I'(s,y(s)) 1 1 ds 

t 

< L S 1 lx(s)- y(s) I I ds 

o ^ 

talcing supremum on both the sides, v/e get 

11 [^U-B]x - [A+B]y I 1 q < L } ix - yl 
which gives a = L. Hence 



a\ 3 
P E 



< 1 . 


Lastly 

t 1 [A+B]x(t) 1 1 < 11 q(t)l L + / 1 !5’Cs,x(s))l 1 ds 

o' 

t 

< a + / ll lx(s)l ds 
0 

and so 

1 1 [A+B]xl 1 ^ < a + ll jx} 1 Q < a + Lr < a +( ^ - |)r = |r 
Thus [A+B] x(t) e whenever x(t) e . 

iind hence by theorem 2.4.5 the given integral equation 
has a unique solution x e 0 with x(t) € . 



GH iPTER ~ THREB 

BOUlDiRY VillJE PROELEI«I 

IICTRODUGTIOE 3.1 

In control ttieory, to find out tiie value 
of the optimal control, hy the help of necessary conditions 
of msjcLmum principle or calculus of variations, one needs 
to solve some non-linear Boundary Value Problems (BVP). 

^so there are some problems which can be directly conver- 
ted into certain form of BVP. Many other problems of 
natural science and engineering quite often lead to B?P. 

So in a natural way this branch of the theory of Diffe- 
rential Equation (DE) has become a very interesting field . 
of research for both the mathematicians and the applied 
scientists. 

Ihere are various methods available in 
the literature for the study of the qualitative properties 
of BVP. Ihe methods [17,18,19 ] which are non-construc- 
tive type, mostly deal with certain class of DE ?/ith 
certain definite form of Boundary Gonditiohs (BG). . Ihey 
are.- based on certain properties of corresponding initial 
value problems, some differential inequalities or some 
assimptions on the uniqueness properties of certain BVP. 
There are few interesting literature [ 2G,21 ] on the 
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non sol'vaoili'ty 1:116 B\rp, iiliese again give resulis 
only for certain class of DE. In the neict section with 
the help of an example [ 22 J we i¥ill show the typical 
benaviour the non-linear BVP exhibits with regard to its 
soliition, 

Jmong the various constructive approaches 
for the study of the properties of the BVP, Quasilineari , 
zation technique [ 10 | is very popular. In this method 
the non-linear BVP is replaced by a. sequence of linear 
problems whose solutions converge to the solution of 
the original problem. Bat's as ::^en.ti one d in the last 
chapter this method does not have anj?" special advantage 
over the methods discussed there, except in the conver- 
gence rate, but at the expense of higher number of 
computations required, Hov/ever we are interested in 
the application of fixed point theorems presented in the 
last chapter for the study •f the existence and unique- 
ness properties of the solutions of the non-linear boun- 
dary value problem (BVP) . 

Somehow there are very few literature 
available in this direction. The booh [ 8 ] gives a 
very good computational approach for the study of the 
most general form of the BVT ir/ith the help of fixed 
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point and other functional analytic theorems. In [ 22 ] 
a certain class of IE with certain type of boundary 
conditions (BG) have teen studied ty GMP. 

The main difficulty in the study of the 
BVP via fixed point theorems is to give a suitahle 
integral equation (lE) representation of the prohlem. 

It seems that the trend in the existing literature 
[ ] is to give the representation with the help 

of Green's function of some linear BTP, In this chapter 
we present a very simple IB representation of the non- 
linear BVP without using the concepts of Green's 
function. •The technique involved in the representation 
is capable tf dealing v/ith v.dde cla,ss of constraints on 
the dependent variable of the PE. Using the theorems 
2*4.5 and 2.4.9 the lEs have been stLidied, Pew numeri- 
cal examples have been presented to demonstrate the 
procedure . 

BOUIIPiPY ViLUE PROBLE&I 3 .2 

In this section we define the BVP vrith 
its various classes and mention some notations and 
definitions. We also give an interesting example to 
show the typical behaviour exhibited by the non-linear 
BVP. Tliroughout this work ?/e shall be concerned vfith 
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ordinary DE. Moreover iiie independent variable V'/ill 
alwa 3 (' s be reo.l. It will usually be represented by t and 
called time. The dependent variable which will u.su.ally 
be represented by z will also always be real and called 
state. Unless otherwise mentioned explicitely any function 
defined over the real line with values in some set will be 
represented by z, whereas x(t) will represent its value 
at t ime t . 

Types of F\rP . 

In theory v/e can define various hinds 
BVP, but in this chapter v/e will be interested only in 
two-point end multipoint cases. 

DSPIUITIOl 3.2.1 

The two point boundary value problem consists of the 
Differential equation 

|| = ±(t) = f(t, x(t)) (3.2.2) 

with the Boundary condition , 

g(x(a), 3 c(b)) = G (3,2.3) 

where f and g are suitable vector valued functions 
and t S [a,b] Cl E-* 
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SIMKFIjB 3.2.4 

I'ollo\’7ing are the two very important , 
special cases of the two point jB7P defined above 

X (t) ^ f(t,x(t), i(t) ) (3.2.2) 

with the boundary conditions 

x(a) = A x(b) = B 

or 

x(a) = A 5;(b) = B 

Tho above two problems have been studied extensively 
in many literature [22] . 

DEOTITION 3.2.5 

The multipoint boundary value problem consists 
of the DE(3.2.2) and the boundary condition given 
by 

h( x( a^ ) j x( ag )y .»..x( ) = 0 (3,2.6) 

where h is some suitable vector function and 
a = a^ < ag. . .< a^_^ < a^ = b t G [a,b] Cl R. 

EXilCPLE 3.2.7 


As an example v/e may have the following third order 
differential equation 



39 


x^(t) = XgCt), XgCt)) i = 1,2,3 

v/itli tiie bourLdary condition given by 

Xi(ai) = 

wliere a = a^< a 2 < a^ = b and t G [a,b] . The above BO 
may talce various combinations of x^ and a^. Bor 
example, _vYe may have 

^ ^1 ^ ~ "*^1 ^2 ^ ^ ~ ^3 ^ ^2 ^ ~ 

DEPIMTIOIT 3. 2'. 8 

By solution of the BVP we Vifill mean an absolutely 
continuous, function 0(t) v/hioh satisfies the BE 

(3.2.2) for almost all t G [a,b] .and the BG 

(3.2.3) or (3.2.6) according to ?fhether the problem 
is tv7o point or the multipoint. 

Throiighout in this chapter, withosit any loss of genera^- 
lity, we will assume that [a,b] =[o,l]CIR and that 
0 “ a^ Sg ... a^ — i . 

Rotations 

will always represent the real 
space of dimension n, will be denoted by R. The 
sjTnbol will denote the class of all continuous 
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fuiLCi:ions defined over [ o . l] i ' B. wiih values in 

0 will be denoted by G. Tlie class of functions in 0^ 
v/liich tables values in S d R^ v/ill be denoted by G^fS"] 

1 .e . 

G^[S] = ( 0 : 0 e , 0(t) G S, SCI S"'^) 

1 1 *1 ln_ will represent the norm in R^, and | .j will 
represent 1 1 • 1 | • We define 

1 lx(t)l 1 = max { lx.(t)[ ; i = 1,2, . . aa} 

i 

where are the components of n-dimensional vector 

33 . 

x(t) . The norm- in G , denoted by [j •! 1 q5 is defined as 

i 1 0 = Sup [ |!^(t)l t e Io,l] 

"fc 

= max t Sup j04(t)jl t 8 [b,!] ” i =1,2, ..n} 

it 

where 0j^(t) is the ith component of 0(t) , 

The norm on the finite dimensional 
operator A(t) in P0'^, a matrix of size n x n, with 
elements a. .(t) 8 0 is sdefined by 

n ^ - 

1 I aI I = max ( 21 Sup j a. .(t ) | ; t b[o, l] , i = lj2, . .n} 

° i 3=1 t 

In a similar way the norm on G(t,s), a matrix of order 
n X n, with elements defined for t,s g[_o,iJ 
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wii^li valu6s in E. is giTBn by 

n 

1 1 &! 1 Q = max { 2 Sup Ig. ;-t,s e[o,lJ 1=1, 2, . , .n ■ } 

i 3=1 t,s ^3 

If "tliere is no confusion possible all the norms except 
the one in R v/ill he represented hy j [ *1 | • 

Ihe symbol will stand for a closed 

44 

sphere of radius r centred on origin^ i.e, 

= (z ; 1 Ixj < r, X e R^ } 


hERIhlll Oh 3^; *9 
Let I and D he open subsets of R and R^ respecti- 
vely and let [ojljdl. Ihe function f(t 5 x) v/ill 
he said to satisfy the property (H) if all t>£ the 
following conditions are satisfied, . 

(i) f(t,x) : I z D R^ is continuous in z e D 
for each fixed t G I and measurahle in t for each 
fixed X G L. 

(ii) There exists an integrahle function m(t) 
such that 

i a.e. in I 

and for all x G L. 

(iii) f satisfies the uniform Lipschitz condition 
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f(t,x) - f(t,y)l L < 


n 




Ml 


foJ? 0,11 31 5 y S D sTid. ci»©, in l, ' 

(iv) f(t,0) = 0 for all t e I. 

liie first tv/o conditions are linovvn as Caratheodor^'" 
conditions and the fourth condition will define a free 
or unforced system. 

3 ., 2. 10 

In this section we present a study on 
the solution propertios of the following BVP feal 


• • 

y + = 0 

y(o) = 0 y(b) = B 


(3.2.11) 

(3.2.12) 


Wo first observe the follov/ing few properties, 
solution of (3.2,11) is a. solution of 

*y - y = O if y(t) < o 


A 


(3.2.13) 


and 


5 r 4- y = 0 if 3 l(t) ^ 0 


(3.2.14) 


The solution of (3.2.13) is y(t) = -a sinlit end of 
(3.2.14) is y(t) - a sint . for some arbitrary constant a. 
low if at t = t^ y(to) = 0 and y'(tQ) < o, then y(t)< o 
iraiiiediately to the right of IJq and hence is a solution 
of (3 .2 .13) , and will remain as solr.tion of (3 .2 .13) for 
all t > t|^ , since -a sinht is always negative . 
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In a similar way if at t = to 7 (^ 0 ) ^ ® then y(t) > 0 
on some interval to the right of t^ and will thus he 
the solution of (3,2,14), This will remain as solution 
upto t = 7 t , at which y(Tc) = 0 with yCit) K o and from 
then onwards it will follow: the equation (3.2.13) in 
a manner mentioned above , 

With this knov/ledge on the solution of (3.2.11) we can 
now analyse the B?r. . 

Case I : h < 71 

In this region we have seen equation (3.2,11) has two 
solutions having both pfsitive and negative values. It 
is to be noted that zero is the trivial solution -•f 
(3,2,13) and (h.2.l4). So in this case we have 

^ a sin t if B > o 

y(t) = s' -a sinht if B < o 

(0 if B = 0 

and hence the solution is unique. 

Case II: b = 7t • 

If B = 0 there are many solutions given by 

{' a sin t for all a > o 
y(t ) = i 

I 0 (trivial case). 

But if B > o there are no solutions. Since solution 
of (3.2.13) always remains negative and of (3,2,14) 
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attains zero at t - tc* But if B < o there exists only 

one solution y(t) = - a sinh t 

where a is to he selected to satisfy y(h) B, 

Case III: h > -re . 


As above if B > o there are no solutions . If B = o 
trivial solution is the only solution, because there are 
no nontrivial solutions which attain zero after t > a . 

If B < 0 , then there will be two solutions. One solution 
is apparent, given by y^(t) = - a sinh t v/ith a constcuat 
a such that y^(b) = B. 

Ihe other solution v.dll be the one which satisfied equation 
(3,2.14) for o t <_ n; and then (3,2,13) for it 1. "t X. b, 
i.e. the one v/hich remained positive in o < t < % , zero 
cit t = 71 , and then negebive for t > 7t. This can be repre- 
sented in the form 




(' c^sin t 


O < t < 71 

■ Co(sinh t - tan % cos3a t) tc < t < b. 

k. ^ 

the constants c^ and c^ should be selected in such a way 
that y 2 (’b) = 3 and y(t) is continuous at t = ti. 

The above simple example shows, how 
complicated the behaviour of the solution of a B\rP may be. 
It also shows that not only due to the large value of the 
Iiip s chit z constant but aZLso because of the length of the 
interval and the BG the EVP may fail to exhibit solution. 
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Tlie applica"tion of OMP 1:0 tlie atove problem will give 
■fcliG roaulis o£ Case I. lliis is because ilie lipscliiis 
condition carries very little information of the function, 
lliis encouraged some people to work on methods not based 
on Lipscliitz conditionj^g2] But so far it has not been 
possible, with the help these methods, to give very 
general results. On the otherhand Ol’EP can. be quite 
satisfactorily Lised to study the very general problem. 
Modifications of norm may also help to give better results. 
But it appears that one will be a.ble to give best possible 
results only for particular cases and not for the general 
problem. 

REPRESEIIlillO h P R O:aLM 3,3 
, We now« give the integral equation repre- 

* 

sentoi;ion the tw» point BVP defined by (3,2,2) and 
(3.2.3) . 

PHE ORffl 3.3.1 

Let the function f(t,x) satisfy the first two 
conditions of the property (H) . Let g(x,y) be 

H 

continuous function mapping D x D into E. . Then 
the Boundary Yalue Problem 

i(t) = f(t,x(t)) g(x(o), x(l)) - o (3,.3.2) 

has the equivalent representation 
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"b 

x(t) = x(o) + g(x(o), x(l)) + / f(s,x(s))ds (5.3.3) 

0 

for oilmost all t e [o,lJ . 

I’JP.9.?' ^*3 are to show that aiiy solution of (3. 3. ,2) 

sOvtisfies (3.3.3). and vice-versa, let 0 he a. solution 
of (3.3.2), then we have 

0(t) = f(t, 0(t)) g(0(o), 0(1)) = 0 

v/hich can he v/ritten a,s 

t 

0(t) = 0(o) + / f(s, 0(s)) ds 
0 


to the right hand side of the above we can add g(0(o),0(l)) 
v/liich is identically equal to zero, without changing any 
thing, thus 

t 


0(t) = 0(o) + g(0(o), 0(1)) + / f(s,0(s))ds 

0 

And hence 0 satisfies (3.3.3), 


Conversely, let 0 he a solution of (3.3,3) 
Then , 

"U 

ri(t) = 0(o) + g(6(o), 0(1)) + / f(s, :6(s))ds 

0 

from above it is clear that 0 is an absolutely continuous 


function and therefore differentiable almost everywhere. 


On differentiation we get 

e(t) = f(t, e(t)) 

also putting t = o, the last IS gives 
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0(<5) = e(o) + g( 0 (o )5 e(i)) 
i.e. g(e(*), 0 ( 1 )) = 0 . 

Thus 0 satisfies equation (3,3,2). q.e .D. 

OORO LLiffi Y 3,3.4 

Let the hypothesis regarding f(t,x) of the theorem 
(3,3.1) he satisfied. Then the multipoint boundary 
value problem 

x(t) = f(ty3c(t) g(x( ) j x( Ug) j . . .x( 3^) ) — 0 (3,3, 5) 

where o = < a. 2 ** - ■^a^ = 1 and g is a, suitable 

vector valued functionj has the eqviivalent represen- 
tation 

t 

x(t) = x(o) + / f (s,x(s))ds + g(x(a^) , . .x(a^)) (3 .3 .6) 

The proof of the above corollary can be given in the 
same wajr aa. has been given for the theorem (3,3.1) . 
Following are few excmples of special cases. 

Consider the B7P in R^ 

i^(t) = f^(t,x^(t), X 2 (t)) 
xgCt) = 3Cg(t)) 

x^ ( o ) ■= a Xg( l) = b 


SXiilPLB 


3.3.7 
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//e can dTorm various functions g(---(o), x(l)) represen~ 
ting tlie boundary conditions mentioned above, Hov/ever 
a straight forv/srd application gives the follovTing IB 
representation 


t 

x^(t) = a + / f^(s 5 X^(s), X 2 (s))ds 


= XgCo) 


t 

t b — Xg(l) + f f £( s jX^ (s) jXg (s ) ) ds 


BXAI\IP1E 0.3.8 

/mother important class consists of the DEs of example 
3.3.7 with the BCs given by 

x^(o) = a x^(l) = b . 

Again we can construct various g functions, but a 

direct application gives 
t 

x.(t) = a+/ f (s,x^(s), XgCs)) ds 
° t 

XgCt) = Xg(o) + b- x^(l) +^f fg(s,x^(s) 5 X 2 (s))ds 

REIABXS 3.4 

The advantage of the IE represent nation given by theorem 
3.3.1 is that it can include many types of conditions 
on the system states, and parameters p) re sent in the 

system, 

Bor example in the case of periodie 

solutions the g function will have the form 
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g(x(o), x(l)) = x(o) - x(l) 

If it; is required "bhcffc "tlie integral of tlie norm of tlie 
solution sliould have some value c then v/e can tahe g as 

1 

g(x) = / 1 ix(s)l I ds ~ c 
0 

The technique of the theorem 3,3,1 can he extended to 
the follov/ing very important class of problems knoim 
£is controllability problem. . 

lEglMTlOIT 3.4,1 

Let I, D, U be open subsets of E, and R^ respec- 

tively with [o,l]CZI . The function f('t,x, u) will 
be said to satisfy the property (P) if 

(i) f: I X L X U ->■ R^ is continuous in (x,u) e D x U 
for each fixed t G I and is measurable in t G I for 
each fixed (x,u) G L x U. 

(ii) there exists an integrable function m(t) 
such that 

j|f(t, x(t), u(t))|i^<.m (t) 

If there exists a function u(x°,x^) = u..G called 
control, such that by using this u the solution of 

x(t) = f(t,x(t), u(t)) 


the LB 
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c£!ii be steered from any given to any given in 
some finite time interval, then we say that the DE is 
controllable. More specifically we use the following 
definition. 

EEEIIIIION 3,4.2 

The controllability problem consists of finding out ■ 

an u e G^[U] such that the BVT 

:x(t) = f(t,x(t), u(t)) (3,4,3) 

x(o) = x(l) = x^ (3.4,4) 

o i 

admits a solution for some given x“ and x in D, 

let g : D X E be a continuo^ls 

fimction siich that 

g (x,y) = ,.0<^x = 0, y = 0 
then the BC (3.4,4) is equivalent to 

g (x(o) - x°, x(l) - x'-^) = .0 
Hence along the lines of theorem 3,3.1 we have 

T HEO REM 3.4,5 

let f(t,x,u) satisfy the property (P). Then the 

Gontrollability Problem dsfinbd by (3,4.3) and (3.4,4) 

has the equivalent representation 

1 

x(t) = x(o)+ I(x(o)-x°,x(l)-x^)+ / f(s,x(s),u(s))ds 

: . ■ . 0 



The proof of the above theorem is very similar to the 
one given for theorem 3,3,1. We avoid repeatation. 

SOLIJlIOh PROEDM 3. '5 

Having given the IS representation we can 
now use the fixed point theorems presented in the 
second chapter to find out the solution of the problem. 

We illustrate the solution problem using theorem 2,4,5 
v/hen the operator B is taken as linear, Ihis will 
simplify the presentation. It is to be remembered that 
B need not be linear. An example with a non-linoajr B 
operator has been included to demonstrate the advantage 
of this type of representation. 

As the BVP (3,3,2) is equivalent to the 
problem (3.3,3), it is sufficient to study the IE. 

Define the operator A by 

t , 

y = Ax y(t) » x(o)-. x(t) +g(x(o) ) ,x(l) )+ /f(s,x(s))ds 

o 

where f(t,x) and g(x,y) satisfy the Jnypothesis of the 

Hr T 

theorem 3.3.1. Ihus the operator A maps C [DJ into 
0^. life are interested to find out a ’f) G C^{ i ] such 

that 

A 0 = G 


v/hich will then mean that 
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■fc 

0 = 0(o) - 0(t) + g(0(o), 0(1)) -i- / f(s,0(s)) ds 

0 

i.e, 0 will Toe the solution of the IE (3.3,3). 


Let us assume a B operator of the folloY/- 


ing form 
y = Bx y(t) 


t 

x(t)- x(o) --Mx(o)- Ix(l)- /Y(s)x(s)ds 

0 


v/here M,E and Y(t) are n x n matrices. It is assumed 
that Y(t) is measurahle a.e . on I and with the norm 
dominated by an integrable function over I. 

With the above choice of B the operator 
[ + B ] talces the follouring form . 


y=[A+B]x 

t 

y(t)= g(x(o),x(l))- !’Ix(o)-Hx(l) + /[ f (s,x( s) )-T(s)x(s)] ds 

0 

Let us define 

L-(x,y) = g( 2 c,y) - Mx - ITy 

E(t,x(t)) = f(t,x(t)) -V(t) x(t) . 

Assume that h(x,y): D x D -*■ satisfies , 

Hh(x,y) - h(u,v)l 1^< Lj^^l |x-uj + Lj^l ly-v| 

and h(0,0) = a. 

also that E(t,x) satisfies all the conditions of the 
property (H) v/itli f replaced by E» 



I^rom tlie explicit expression for tlie 
opero-tor [a + b] we see ttiat for any x G [ D ] ,y 
is an absolutely continuous function. In a similar vmy 
for any given absolutely continuous (AG) function y the 
operator B shows that x is an AG function. Bote that 
in B operator y will always be given as an AG function, 
because according to the theorem 2.4,5 this will be 
generated by the operator [ A + B ] . 

Differentiating the expression for Bj 
which is possible because of absolute continuity of y 
and X, we get 

x(t) = V(t) x(t) + y(t) 
while putting t = o, we find 

Mx(o) + hx(l) = ~y(o) 

How if the matrix [ M + H ^(l,o) ] is non singular 
where ^tjt^) is the fundamental matrix of the system 

x(t) = Y(t) x(t) 

with = I > 0 - t^ < t < 1 

then by Lemma 1.5.6 , the inverse of the operator B can 
be given by 


1 
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where the Green’s matrices H(t) and G(t,s), as mentioned 
in section 1,5, are given hy 

. “1 

H(t) = $(t,o) [ M + I ] 

G(t, s) = H(t)M ^ (o, s) 0 1. s < t 

= ~H(t) I® (l,s) t < s < 1 

To compute the value of p for the operator B we rewrite 
the inverse ejq)ression as 

x(t) = - Ii(t)y(o) + .G(t,l)y(l) - G(t,o)y(o) + 

+ /^&o('*^js)y(s)ds 

° 1 
= - H(t)[ M+I 1 y(o) - H(t) B y(l) + / G„(t,s)y(s)ds 

0 

Where 

Assume that 

with the above hounds we compute that 

I l^llo i ^2+ 

and hence 

P = P'1 + 
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To calculate tlie constant a fon tlie operator [ A+B j we 
observe that, for any x,y 8 C^[ I) | , 

11[ A+B] x(t)^ [ A+B ly(t)ll^ 

< j lli(x(o),x(l))-}i(y(o),y(l))| +/ j [bCsjxCs) )-I'(s,y(s) ) | j^ds 

o 

t 

lx(o)-y(o)| 1:^+I^1 jx(l)-y(l)| 1 ^+ 1 ;^ / 1 1 x(s)-y(s) ] 1 ^ds 
and lien.ee 

11[ i+B ]x - [ A+B jyll^ 1 +L^)\\^ ~ yllc 

therefore a can be taken as 

a = + Lj, 

Let be a closed subset of D. Then for any x G ] 

v/e have 

11 [ ^B]x(t)|l^ 

1 

<1 lh(x(o),x(l))-b(o,o)l Ij^+I lli(o,o)| 1^+/ Lj.1 lx(3)l|j^ds 

< lx(o)l Iji + Ijigl l^a)l Iji + 1 l^-l In ln'^3 

or • ■ 

1 1 [ i-+B 1 Q 1 I I c I 1 ^1 1 c 

< ar + 1 1 aj 1 ^ • 

Tims 1}o satisfy coiidition 

[ iH-B ] Sj,/p . 



v;e need 
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+ 1 i I c i 
1 I ^1 !c i - ap)/P 

Hence we liawe proved dlie following "blieoreni. 

lliEORM . 3,S.l 

Suppose the following conditions hold 

(i) The function ^(tjx) defined hy 

^(tjxCt)) = f(t, x(t)) - V(t) x(t) 

satisfies £ill the conditions of the property(H) 
mentioned in definition 3,2.9 with f replaced by H. 

(ii) There exist nxn matrices Li snd H such that 

the function h s D x D defined by 

h(x,y) = g(x,y) - llx - Hy 

is continuous and satisfies 

I |h(x,y)~ h(u,v)| 1^ < lx>-u| jy-vl 1^ 

ifiii) There exists nxn matrix Y(t) vdth integrable 
elements such that the matrix [ M+ H$'^ (l,o) ] Is non 
singular y where the funaamental matrix 

of the system 


i(t) = v(t) x(t) 
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with Mto,to) = 1 , 0 = < t < 1 

(iv) The functions H(t) and G(t,s) as defined below 

H(t) = ?>(t,o)[ M+ IT^(l,o) ]“^ 

G-(t,s)= H(t) M ^(o,s) 0 s < t 

= ~H(t) E'^’(l,s) 1 

satisfi^ the hounds, given hy 

nH(t)[ Ml ]| l^< 1 1HC1;)U||3 < Pg llSg(-t,s)ll^< Pj 

vjbere (}g(t,s)= j| 6(1;, s), suoli that 

< 1 

where a = and p = p^+ 

(v) The element h(o,o) = a has the hound given hy 

i ta|lc i ^(1 - ccPVp 
for some r > o such that S^cv 1. 

Then the 5VP 

x= f(t,xCt)) g(x(o), x(l)) =0 

has a unique solution ^ G 0^[ ] . 

ElLaiPLES v3.6 

In this section we present few numerical examples on 
the non-linear BVP. The idea is to show that the 
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..-(.aijUiaptlons on i:heorera 2,4,5. and "blie reprerien'fca'tion 

oy tlieorem 3,3,1 are valid and tlae set of proTolems 
tlmt 0821 be solved is nonempty. 


EXIMM 3.5,1 

It is required to find out the solution of the following 
problem 


(t) = 


¥ 


i7xf(tT 


1 ^i( ^) 

^ 1+ x|(t") 


Xg (t) 

x^(o) = 1 XgCl) = 1.5 


Solution ; the theorem 3.3.1 the above problem 

has the equivalent representation given by 

"fc 

(t) = 1 + / I 2 Cg(s) / (1+ X 2 (s)) ds 

Xg(t) = 2Cg(o)+ 1.5 - X 2 (l)+ / ^ X^(s)/(1+X^(s))ds 

Ihus the operator A will have the form 

■fc ' 

y = liX TqC^) = l'-2C^(''^) + 4 Xg(s)/(1+Xg(s))ds 


(t) = Xp(o)+1.5 -XgCD-XgCt)^- I /x^(s)/(l+xj(s))ds 


^2 


natural choice of B is 
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y = Bx y^(t) = x^(t) 

y2(t) = Xg(t) - XgCo) + ^2^^^ 

Thus [ -I + B ] v/ill have the form 

y = [«-.+ B]x y^(t) = 1 + i / X 2 (s)/(l+ x|(s)) ds 

72 (t) = 1.5 + i / x;^(s)/(l+x^(s))ds 


Both B Sind + 3 ] are defined over the \vhole space 

P t* 1 

C , .'md hence the condition [ ^+B ] will 

not be necessary- 

The operator B has the inverse given 


-1 

= % : 


= y^Ct) 

XgC'b) = y2('t^) + “ ^2^^^ 


Then for any y G 0 , we have 


x. (t ) 1 


i 


XI 1 = max Sup I lx(t) I j 

° i t 


n 


max Sup I 


< max Sup 
“■ i t 


_ _ i ^ ° I 


max Sup 

3ly2(-t)l| ^ ^ 


0 


3 I 


j|x2(t)j 

ly2('fc)t 


< I 


1 0 


.0 3j 
and hence ^ = 3 . 


• Hylic i " Ihll 
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ap = 1 .3 = 3/4 < 1 

the above problem has a unique soMtion. 

We have the follov/ing numerical values ior x^(t) and 
XgCt) as the solution of the BVP; . 

t = 0.0 0.2 0.4 0.6 0.8 1.0 

X ( t )= 1.000 1.024 1.047 1.071 1.094 1.117 

X ( t )= 1.375 1.400 1.425 1.450 1.475 1,500 

2 
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e:;oipie 


,3.6.2 


Cn,. -.iic.r :iOv/ tiio follov/ing oxaiaplo / 


(t) 




1 XgCt) 

‘i l+x|(t) 


X^(o) = 1 


X.,(t) = Xp(t) + 


1 X^(t) 

4 l+^(t) 


Xp(l) = 1.5 


3ol^.ion : Define the operators in the following way 

y,(t) = l-x.,(t) - /x (s)ds + ^ f X2(s)/(l+x|(s))ds 


B 


t 


>^2 


p(t) = 1.5-Xp(l)+ XgCo)- X2(t)+ / XgCs) ds 




+ J / x^(s) / (l+x^(s)) ds 
t 

= x,(t) + / X.(s) ds 

0 

t 

=-. x„(t)- XgCo) + X 2 (l) - / -g(s) d£ 


a-3 y^Ct) = 1+ I /“xgCsVCi + x|(s))ds 

y,,(t) = 1.5^- ~ / X^(s)/(1+ 


0 ' 


Cl,. 


arly all tto operatorp are defined and oentimoup o^er 
compute the value of ? and cc and show 


v/liolo ' of , E- ' • „ 
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t licit "^1) tlien ttieorem 2.4.5 v/'ill assure the existence 

and uniqueness of solution. 

I'roM the first equation of the operator B we get 


x.(t) = -X. (t) + y. (t) + , X. (o) 


Ji (o) 


The solution of which is given by 


X. (t) = e”V-,(o) + / e“ y. (s)ds 


y.('t) - / e‘ 
0 


y^(s)ds 


1^1 (■t)l < Sup ly.(t)l + / e" ^®|y^(s) |ds 


< Sup [ 1+ / e 
" t‘ 0 


] . jy^l < 2iy^| 


Similarly the second equation is equivalent to 
XgCt) = X2(t) + ygCt) xgd) = VgCo) 

the solution of which is given by 

Xg (t) = e^"^ ygCo) + / y2(s)ds 


'[ ygCo)- 73(1) ] + ^2^’^^'^ { 
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The above gives 

l^sio = sup iXg(t)| < Sup [ 8e'^-ii+ /e^-s del ly.l 

u ij ^ 

= Sup [ 2+ j.jygl^ < Siy^j^ 

Prom "tiio alDovo "two ineqi.iali'ties we gei} 

ll-^lloi 3||yll^ 

and hence p = 3 . 

Proceeding along the lines of the problem 3.5.1, we 
compute the . quantity a for [ n+B ] as 1/4, giving us 

ap = 4 • 3 = 3/4 < 1 

and hence the BVT has a unique solution. 

The following table gives the numerical values of the 
solution ■ 

t = 0.0 0.2 _0.4 0.6 0.8 1.0 

x^(t) = 1.000 0.8372 0.7064 0.6007 0.5143 0,4428 

Xg(t) = 0.4548 0.5922 0.7589 0.9607 1.2050 1.5000 

EOlPLE 3.6,3 

We solve the following problem with nonlincax integral 
B operator using the theorem 2.4,9. 
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yir 


4+ x: 


- X, 


x^(o) = 0 


X, 




- Sx. Zp(l) = -1 


Solu'tion I In tliis case -the operator A is defined as 


p u 

/;], ("b) = + f ^2 

yg(t) = -XgCt) 4- XgCo) + 1+ XgCl) 4 

"t/ "ij> 

4 / X (s)/ (44 X ^(s))ds - / 2x^(s)ds 

lake the following non-linear B operator 


U 

B yi('b) = 2t^(s) h£ 


2 

3rg(t) = XgCt) - XgCo) - Xg(l) 4 / 2x^(s)ds 

which gives [ A 4 B ] as 

t 2 

^U-B y.(t) = / XgCs)/ (4 4 X2(s))ds 

o 

t 2 

^^(t) = 1 4 / X^(s)/ (4 4 x^(s)) ds 

Before finding out the values of the parameters a aud ^ 
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let us calculate the Lipschitz constant and the hound 
of tlx. function x/ (a+x^). 

(d/dx) (x/(a+x^)) = (a-x^)/(a+x^)^ 

v/hich attains the maximum value at x =0. Hence the lip 
Lipscliitz constant is l/a. The function x/(a + x^) 
attains its maximum at x = fa and thenofore the hound 
is fa/2a = l/Zfa . 

?/q find out the inverse of B. from the firrit equation 
of the operator B, it follov/s that 

2:^(0) = 3^1(0) = 0 

v/hich has the solution 

x^Ct) = / y.(s)ds 

y^(o) = 0 follows from the first e nation of [ i^-B ] . 
Similarly the second eqnation gives 

yg(t) - 2x'J(t) = yg(t)- E[ /e-^'^-=>yj^(s)<ls ] 

XgCl) = 72^°) = ^ 


from which we get 
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t s / ^ ^ . 2 


XgCt) = 1+ 72^^^ ■“ e'*'^®”'^^ y^(v7)dw ] 


ds 


Hov/ subs'fci'tu'fce "the values of 7^(1^ ) and y^Cf) fncm tlie 
operabor m-B ^ in tlie above two equations for x^Ct) 
and XgCt) and get the following set of lEs. 


x.Ct) = 

^0 4 +X2^(s) 


XgCt) » 1+ / --ig- ds - / ™~V 
o 4+ x^ 0 4+xJ 


2 /{ 


1 o 


4+ X, 




_ dw y ds 


1 - / 


1 X. 


t s / \ JS-o ^ 

ds - 2 / [ / dw] ds 

1 0 4 




t 4 + x' 


Ihe solutions of the above two e qu at ions can now be 
studied by GICP. Since they are defined over whole of 
we need only to test the first condition of theorem 

g 2 

2,2.1 . het us define the operator 1 : 0 - C by the 


following 


T 


r 1 

X 


x^ 


(t) 


/ X2(s)/(4 + 


ij s ^ 

- A,/(4+ 4)1= - S/[ f Ml® 


^4 



Which gives 


j [Tu - Ivj = max Sup j [t 

t 


■u. 

\y 

1 

(t) - 1 j^lj 


(^)ll 


< max Sup 
t 


It 


J 1 1^2(s)/(i-+U2(s))- Vg(s)/(4+Vg(s)) | [ [hs 


1 ■^-1 V 1 S U^ Vo . 

/| — 2 ds + 2/j j / — K- + — ““9 I I 

o 4+u^ 4+v o 0 4+Up ^+”^2 


Up Yp 

I j dwds 

4+ uj 4 + Vg 



0 1/4 

1 

|u^(s) - V^(s)[ 

< max Sup 1 1 


!l • / 


- t 

1/4 1/4 

0 


ds 


< t 11^ - ^l.c 

Hence a = l/2 < 1 which, hy CaiP, shows that the 

opercitor T has a unique fixed point. In this problem 
the sequence {x^} is generated as follows: 

4(s)/(4+(4(s))t ds 

X r\ 


n 


t ji 2 


x"^(t) = 1- / ds - 8 / [ x“(s) ] as 

* 4+(Xj^(3))^ ^ 
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st.-.rting v/ith any arbitrary x°(t) aaid x°(t), 

The table below gives the numerical values of the 
solution functions x^(t) and Xg(t) 

^ = 0*0 0.2 0.4 0.6 0.8 1.0 

Xi(t)=0.0000 -0.0359 -0.0654 -0.0896 -0.1094 -0,1257 

X2(t)=-»9954 -0.9960 -0.9975 -0.9988 -0.9997 -1.0000 

3.7 

It has been observed in the last section 

that the choice of the operator B has been dictated by 

the form of the operator A. In many cases it may not 

be so. It will usually be simpler to select a lineae: 

B operator of the form given in section 3,5. The choice 

ff tl^ matrices and V(t) in the theorem 3.5,1 should 
0 

be made in such a way that the quantities ^h2’^B 

are reduced to very small values. Bor example the 
matrix Y(t) ma37- be selected in such a way that the 
Brechet Derivative of the opera,tor 

B(tjz) = f(t,x) - V(t)x 

over a certain region inside D is minimised which then 
v/ill ensure the lovi/est value of Dj, . 

In example 3,6.1, because of the very/lovr 
value of the Lipschits constant of the function f (t ,x) 



it beoii possible io select an algebra-ic 


B operator 


Tliiu is an advantage of tlie representation presented 
lit/re over tiie nietiiod discussed in 8 ^ wiiicli will 


r.lwr.ys give an integral B operator and will Be dictated 


by tlio esjpression of the operator A. 


It is interesting to observe in the above 
examples, that, if the eleraent y, in the range of the 
operator [ ] is substituted c in the domain of B“*^, 

the resultant expression will coincide ?/ith IE represen- 
tation of the BVT as has been presented in ® 
course tiiis will be so only when B is linear. 


The essential idea behind the fixed point 
theorem 2.4,5 or 2.4.9 is to extract- a portion of the 
operator A, v\fiiich prevents the direct application of CMP, 
and find its inverse in;,- terms of the remaining part of 
A and then apply CEIP. Examples given in the last 
section, specially the example 3. .6.3 illustrates this 
procedure . 

A little observation of the examples v^-ill 
show that the direct application of Gl'EP to the operator 

A+ B] will usually give less restrictive conditions 
on A then will bo roq^^ired by setting ap < 1 where 
^ a = lips chit z constant of [ A + B ] . This 

conclusion is the converse of the one .observed in 8 ] . 



OHiPTER - POUR 

GOIITROLIABILITY PROELM 

I ITRODUGEIOE 4.1 

Pormally spealcing a system z = f(x,u) 
defined over the real space, v/here f,x,u are suitable 
vectors, is said to be controllable, if there exists a 
control function u , \?hich steers the system from any 
given 2c(tQ) = Xq to any given x(t^) = Mer since 

the notion of controllability was introduced by Kallman, 
various authors have v/orhed on this problem in various 
different waj/'s and a lot of materials are available in 
this direction, Pollowing is a brief discussion of some 
of the interesting literature. 

Controllability for linear systems is 
well developed, any standard teirb book on control theory 
may be consulted for this result. In [ 23 ] a good 
chapter has been presented for the study of non-linear 
controllability problems. In one method, it has been 
shown that the domain of null controllability, i.e, the 
set of points that can be steered to the origin, is an 
open set under the assumption that the system when 
linearised at origin is controllable. Another result 
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v/liich gives global coirtrollabili'ty , malces us© of flie 
Lyapunov function and 111 © confrolla'bili'by of fiie sysfeui 
linearized af fhe origin. In both, of the above results 
it has been assiiuied that the control is restricted in 
some subset of the -.underlying space. 

Some results are available regarding the 
use of fixed point theorems for the study ©f contro- 
llability properties, lu [ 24 I Schauder’s fixed point 
theorem has been used for the study of ■ 

• 

X = A(t,x) X + B(t,x) u . 

by representing it in the form of a suitable integral 
equation, which makes use of the corresponding linear 
system x = A(t,z)x + B(t,z)u for some z belonging 
to certain class of functions. Global result has been 
obtained under the assumption 

|a^.(t,x)l < M lb^^(t,x)[ ■< 1 

for all t, X, where a^^ and b^^ are elements of the 
matrices A(t,x) and B(t 5 x) respectively. 

The controllability property of 

X = g(t,x) + h(t,u) has been investigated in an 
interesting peiper [ 25 ] with the help of fixed point 
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theorem of Ky Pan [ 26 ] . With some convexity ass-unip- 
tions on the domain of i(t) controllability results 
were ootained, while without this assumption an 6— 
controllability has been derived which assures that the 
final stake will be within any given G— sphere around 
the required given state. 

Ihe above methods are essentially 
qualitative methods, in the sense that they do not give 
any procedure to compute a control that will be required 
to steer the system. In [_ 27 ] , a constructive method 
has been presented for finding out the steering control 

m 

of the system x = 9(x) + Q(t) where Q(t) is a vector 
v.dtli only one non-zero element. Ihe method generates 
a sequence of controls by successive use of the equi- 
valent integral equation representation of the linear 
system x = A(t)x + f(t) and non-linear system 

= e(x) - A(t)x + Q(t), which converges to the required 
control, ihe control is selected from the class of 
piecev/ise constant functions over a. given control interval. 

In this chapter, we present a study of 
the controllability property of the system x = f(t,x,u) 
by the help of the fixed point theorem presented in 

chapter two. The controllability problem which is a 



73 


Icijid of bouii'iarj'" value prolleiD. lias leen represented 
by an equivalent integral equation and tben investi- 
g&.ted« vAn e^csriple lias been included to illustrate tlie 
procedure. 

4.2 

Let X and y be t'wo vectors, tiien 

[ xy ] will denote a vector each, component of V7hich 

is the product of the corresponding components of the 

2 

vectors x and y, [xx. ] vail be denoted by [x] 

We follow the notations of section 3.2 regarding the 
norms etc. 

DEglhlllOh 4.2.1 

The function f will be said to satisfy the 
property (PP) if the following condition hold 
(1) f(t,x,u) ; I X D X TJ - where I, D, U 
are open subsets of R, R^ and R^ respectively, 

■Titli [ 0,1 ] dJ • 

(ii) f(t,x,u) is continuous in (x,u) G L x U for 
ea.ch fixed t G I and is measuraole in t G I for 
each fixed (x,u) G D x U. 

(iii) There exists a Lebesgue integrable function 
Bi(t) such that 
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1 I x(t,x,u) j <„ ni(t) a.e. ini ajid for all 
(xj-u) e D X TJ. 

(iv) f saiisfios "blie unifono. Iiipsch.i'tz condition 
in D X U, 

jjf(t,x,u) - 1 < L^l jx-yj l^+lgl |u - vj 

for almost all t 8 I. 

(tt) f(t,0,0) = 0 a,e . in I. 

Througliout in this chapter we will consider the follow- 
ing dynamical system 

x(t) = f(t,x(t), u(t)) (4.2.E) 

DEP IITIlI Oh 4.2,3 

A hounded measurable function u with values in 
U v/ill be called an admissible control. 33^ traje- 
ctory we v^ill mean any absolutely contimioxis (AC) 
function 0 that satisfies the system 4,2.2. with 
some admissible control. 6(t) v\dll denote the state 
of the system at time t. 

pE mill Oh .4.2.4 

The system (4,2.2) is said to be controllable if 
for any given set { [ ] Cl 1 , a, ,b 8 E ^5 there 

exists an admissible control such that the 
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tra.leo-fcor^? of (4.2.2) satisfies the condition 

^ = "b (4.2.5) 

i’or the s alee of simplicity and without Ipss of any 
generality we v/ill assume that [ t^ ] = [ o , 1 ] . 

Prom now onwards hy controllahility problem we ?fill 
mean the follo¥/ing Bounda.rj’- Talue Problem 

x(t) = f(t,x(t), u(t)) (4.2,6) 

x(o) = a x(l) = b (4.2.7) 

For some a,b belonging to D • 

RmBSEI llATIO h 4.5 
THB OBEM 4.3.1 

Let f(t,x,u) satisfy the first three condition of 
the property (PP) . Let K be an m x n matrix 
such that !<;[ x ]^ = 0 ==> x = 0 for any vector 
X G D. Then eny solution of the following IBs 

t 

x(t) = a + / f(s5x(s),u(3))ds (4.3.2) 

0 

t 2 

u(t) = u(l)+ K[ a^b+ / f(s,x(s),u(s))ds ] (4.3.3) 

0 

is a solution of the controllability problem (4,2.6) 
and (4.2.7). 
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Pr(30f : Let 

satisi^r the system of 
Then (4,3.2) gives 


V e u j end e e o^[ d ] 
equations (4.3.2) and (4.3.3). 


t 

G(t) = a + / f(s, e(s), v(s)) d-s (4.3.4) 

0 

Clearly 0 is an AC function, hence hy differentiating, 
we get 

0 (t) = f(t, e(t), v(t)). 

and at t = 0 , (4.3.4), gives 

0 (o) = a. 


Setting t = 1, in (4.3.3), vre get 


/ , 2 

v(l) = v(l) 4 - K[ a-h + / f(s,e(s), v(s))ds ] 

o 

1 2 
i.e. K[ a-D + / f(s, 0 (s), v(s))ds 1=0 
0 

which, hy hypothesis on K, gives 
1 

a + / f(s,e(s), v(s))ds = h . 

0 

But because of (4.3,4) we see that the left hand side 
of the above is equal to 0 ( 1 ), and hence 


0 ( 1 ) = b 

Thus 0 and V satisfy (4.2.6) and (4.2.7). 


Q » S » X) 



77 


SOLUTIOl PRQBLM'I 4.4 

Tiieorem 4,3.1 shov/s ■feh.ai; ins'fcead of studying tlie 
controllability problem defined 13^ (4.2.6) and (4.2.7) 

VYe can study the lEs (4.3.2) and (4.3.3). let us define 
the operators A and B by 

t 

^ + / f (s,x(s) ,u(s))ds 

o 

t p 

3rp(t) = u(l) - u(t) + Iv[ a-b+ / f(s,z(s),u(s))ds ] 

o 

t t 

B y^(t) = x(t) - /V.(s) x(s)ds - /W. (s) u(s)ds 

t t 

ygCt) = u(t)-u(l)~ K / V2(s)x(s)ds - K / ¥fg(s)u(B)ds 

0 0 


Hence the operator [ A+B ] will have the form 


AfB y. (t) = a + / [ f(s,x,u) - V X - ¥^u ]ds 

o 

. t 2 t t 

yp(t) = E { [a^b + / f(s,z,u)dB ] -/ VgZ ds -/ W^uds 

t Q o 0 


With the operators defined as above, we can proceed 
along the lines of section 3.5 to conclude a theorem 
like the one presented there . Boilowing theorem is 
presented for a simplified case when x and u have the 
same dimeiisio]::^ . 
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Let us denote by and 0^ respectively the identity 
and the null matrices of size nxn. 

THEOR EM 4.4.1 

Let f(t,x,u) be a mapping of into R^ and 

a,b be any two elements in R^. Suppose that the 
follov'/ing conditions are satisfied 
(i) There exists two nxn matrices ‘V‘(t) and W(t) 
vdth the norms bounded by integrable functions 
such that the matrix [ M + H <5(l,o) ] is non- 
singular. Y/here 



a 2 id o(t,tQ)is the fmidainental matrix of the system 



Y(t) 

\?(t) 

z(t) = P(t) z(t) 

P(t) = , . 

Y(t) 

w(t) 


vdth <l:(tQs'to) = ^2n ’ "^o = ° * 

(ii) The function E defined by 

E(t,x,u) = f(t,x,u) -.Y(t)x - ¥(t)u 

satisfies the property (PP) with f replaced by P 
and I,D,U by R, R^, respectively. 
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(iii) The Green’s inatrix G(t,a) defined by 
G(t,s) = 5(t,o)[ M+II Kljo) ]“^M §(o,s) o<s<t 

= <^(t,o)[ (l,o) s(i,s) t<s<l 

has the norm 6 , i.e, j | Gj ] ^ = 6 , such that 

aCL^ + Lg) <1 

Then the controllability problem defined by 

x(t) = f (t ,x(t) ,u(t) ) x(o) = a , x(l) = b 
has a solution. 

Proof : Suppose that the assumptions hold. 

Py theorem 4.3 .1, the above controllability problem can 
be solved using the lEs 

t 

x(t) = a + / f(s,x(3), u(s))ds 
o 

t 

u(t) = u(l) + a~b + / f(s,x(s), u(s))ds . 

0 

Accordingly let us define the operator A by 

t 

A y.(t) => a-x(t) + / f(s,x(s),u(s))ds 

^ 0 

t 

yo(’b) = 11 ( 1 ;) + f f (s,x(s) ^ 11 ( 3 ) ) ds 

^ ' O' ■ 

and select the operator B as 
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t 

B y^Ct) = x(t) - / VCs) x(s) is - J W(s) u(s) da 

^ 0 

/ \ “t 

ygCt) = u(t)- u(l)- / V(s) x(s)ds - / W(s) u(s) as 

0 0 

Then "the oper-ator [ A+B ] talces "the form 

t 

-d -3 y^(t) = a + / B(s,x( 3 ), u(s))ds 

0 

t 

y 2 (t) = a-h + / P(s,x(s) ,u(s))ds 
0 

It is easy to see that the operator B is equivalent to 
the following DE 


x(t) 

■ 

v(t) 

■y(t) 


x(t) 


y^C-t) 

u(t) 

1 

1 

T(t) 

w(t) 

1 

u(t) 


y2('fc) 


with the BO 


On 


x(o) 


On 

- 

On 


x(i)* 



r 

o 

H 




u(o) 

+ 

On 

-In 


u(l) 


ygCo) 


The solution ff which is given by 


—> — 1 
x(t) 

; = H(t) 

'yi(o> 

1 

+ / SC^jS) 

ygs) 

u(t) 

«• J 

y2(o) 

0 

i 

y2(s) 

L 


— 1 

v/here H(t) = ^(t,o) [ M + H $(ljo) ] 
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Substitutiiig the value of y^(t) and yg(t) from the 
expressions for the operator [ A+B ] in the above IS, 
we get a nev/ operator I . given by 


x(t) 


a 

1 

F(s,x(s) ,u(s)) 

u(t) 

= H(t) 

j Or— t) 

+ / G(t,s) 

0 

S(s,x(s),u(s)) 

L. ^ 


- 


L 


(4.4. E) 

Clearly T is a mapping from 0 ^ into 0'“' satisfying 

the Lipschitz condition given by 


l|Tp - Iillc i «•( I'l + H >lli’ - allc 


• E n 

vfhere p,q. 6 0 ' . Since ^(I^+Lg) < 1, by GM£ I has a 

fixed point. Therefore by theorem 2,4.9, because T is 
nothing but B~ [ A + B ^ there exists (x , u ) such that 



* X) . 

WIW S 4.4.3 


(i) As has been mentioned before, the advantage of the 
method used here is that it will help us to find out the 
control function numerically. The method presented 
here finds out only one unique control, whereas in 
reality there may be many controls that can be used to 
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steer the system. This has happened because the IE3 
used, to find out the control, puts an equality constraint 
on the control vaniahle, 

(ii) The idea of augmenting the IE (4.3.2) by (4.5.3) 

is to get a composite operator mapping into itself, 

thus removing the difficulty in handling IB (4.3,2) ?;hich 

maps' u into C • • 

(iii) The IE (4.3.3) may be replaced by its most general 

form 0 = h(u(t), t) + g(x(t), b) where h: U x I 

in. n 

g : D X D -* E and 0 is the null element of E , with 

the assumptions that 

h( u(l),l) = 0 and 

g(x(l), b) = 0 x(l) = b . 

Hovrever it appears that the particular form presented 
in theorem 4.3,1 is most suitable for general problems, 

EIIIIPIB 4.5 

We '^resent the following example illustrating the 
applicsition of theorem 2.4,5 in controllability problem. 


x(t) = G 


x(t) 

1+ x2(t) 


_ - u(t) 

1+ u^(t) 


x(o) = a 


x(l) = 0 . 
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Solution : Theorem 4.3.1 gives the follovfing IB 

roproGentation of the above problem. 

t X t u t 

-(t) = a + / e ds + / e — — - ds - / u(s)ds 

0 1 + Z*^ 0 1 + 0 


. t X t U’ t 

ii(t) = u(l)+ a+ / e — + / e — ^ ds - / u(s)ds 

0 1+z'^ 0 l+n* o 


Hence v/e con define the operator A a.s 


y-iCt) = a — ■ 3 ^(t) + S (G . —p-.. + G — n) ds 

^ 0 1 +:^ 1 + u* 


yp(t) = u(l)- u(t) + a + / (6 + G g - u)ds 

^ 0 1+ 1+u 


Select the following B operator 

t 

B y 2 ^(t) = x(t) + / u(s) ds 

t 

yg(t) = u(t) - u(l) + / u(s)ds 
Thus the operator [ A + B ] tehees the form 


44 .+ B 




t 

=: cl 4- G 

X 

- --“p + 

G 

U 

c 

0 

1+x^ 


.1+u' 

+ 

ll 


G 

u 

0 

1+zir 


1+U' 


ds 

ds 


To compute the inverse of 


B we observe that B is 
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n;j_iiivalent to 



The fundamental matrix corresponding to the matrix 



is 






The Green’s matrices are given hy 

-1 

H(t) = ^(t,o) [ M + IT $(l,o) ] 

1 e^"'‘^(l-e'*^)/(e-l) 

0 el-V(e-l) 

G(t,s)=: J G^(t,s) = H(t)M$(o,s) 0 <_ s < t 


where 

G^(t,s) 


1 (e"-l) + el--^+".(l-e'‘)/(e-l) 

0 
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0 (l-e^/(e-l) 

0 e®-V(e-l) 


'Tliws the solution is given by 


x(t) 

= H(t) 

T-,(o)' 

1 

yps) 

u(t) 

1 

+ / G(t,s) 


72(0) 

0 

yg(s)_ 


ds 


G(t,l) 

XAi) 

1 

y-, (s) 




X 

- / 

1 

ds 




■ ygd)^ 

0 

ygCs) 



ygCt) 


where we have used the fact that Gr(t,o) = H(t) 
It can be easily shovai that 


I 1 < e : 1 |Gg(t,s)| [ < e 

for all 't 6 [ 0,1 ] for all t,s G [ c,l ] 

With the above two estiiaatos vre find that the Viilue of 
^ = 2e 4- 1. Irom the expressions of [ + B ] ?/e can 

ea.sily find that its Lips chit z constant a = 2G . Since 
[a + B ] is defined over the ^vhole space, it can be 
concluded that the system can be steered from any 
x(o) = a e E to the origin in the thne interval [ o,l ] 
provided G is such that 


EG (2e + 1) < 1 
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The table below gives the result of a sample calculation 


for C = 0.02 and a = l 




t =0.0 

0.2,:-.-. 0.4 

0 . 6 

0.8 

1.0 

x(t) = 1.0000 

0.7187 0.4363 

0.2935 

0.1333 

0.0000 

u(t) = 1.5670 

1.2900 1.0640 

0.8780 

0.7244 

0.5974 

ikn interesting s 

imp lor result can be obt 

a.ined if 

we 


select the IE representation of the following form 


, V t t t 

x(t) = a + / e ds + / e ds - / u(s)ds 

0 0 l+u*^ 0 


t , -n . 

0 = (l-t)a + a + / e — '—p ds + / e •" — ^p ds - / u(s)ds. 

o l+x^ 0 0 


Define the operators in the follovj-ing way 


^ y-t('t) = a-x(t) + / e — + e —-p - ds 

0 l+x'" 1+u 


y- (t) = (l-t)a + a + / 8 — p + e — - u ds 

^ o 1 + X l+u2 


t 

B y.i('fc) = x(t) + / u(s) ds 
^ 0 

t 

yp(t) = (l-t)a + / u(s) ds 

0 
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y (t) = a + / e + e 

o l+x^ 



x(t) = (l-t)a 

The equation for u(t) \ 7 ill have solution if G < 1. 

Note that the trajectory’' is a straight line joining 
the point a 'vith the origin. 

lollov/irg tahle gives the control for a= l/2 and G = l/4. 

-fc = 0.0 0.2 0.4 0.6 0.8 1.0 

u(t) = 0.7185 0.7039 0.6854 0.6632 0.6381 0.6112 
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CliiPTER _ 

STAEIIIirY PROHLEBI 

ISYROMCTKi; 5.1 

In tills chapter we study the lollov/ing dynamical system 

X (t) = f(t, x(t), u(t)) 

tijc v i. n iahls s j t^x^u^ always rcalj 'will "be ca.lled as 
time, state and control respectiveljr . We call the control 
u, as open loop control when u is dependent only on t, 
and as close loop control when it is a function of both t 
and X. 

Ye present the stability properties of 
the above system with closed loop control using Banach’s 
Eixed point theorem, and with open loop control using 
a modification [ 29 ] of Kalcutani’s fixed point theorem[3] . 

Eixed point theorems of Banach, Schauder 
have been used [ 3l ] for the investigation of stability 
properties. A good number of interesting rersults ma,y 
be •btoined in [ 28 ] , In ['30 ] Ij^choiioff ’ s theorem 
has been used with the help of some differential 
inequalities. In all these cases the system considered 
v/c.s unforced i.e, vltliout the control function. We 
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atudy a more general problem in some what, different way. 

i 

So far as the knowledge if the author 
goes, theorems, liJce Kakut ani * s , have not been used for 
the solution of stability problems. However Ky i’an's 
theorem [ 26 ] , which is a genersHization of Kakutani’s 
theorem in locally convex linear topological spaces has 
been used [ 25 ] for the investigation of the controllar- 
bility problems. One difference we have here in our 
method is that we dcnot require any convexity condition 
on the system, ihe method presented here can easily 
be, without very little change, extended to the case of 
closed loop control. 


Consider the DE 


HEEI IITI OH 5.2 


x(t) == fCt, x(t), u(t)) (5.,2..l) 

with f(t,0,0) = 0, 

defined over the real spaces, where f,x,.u are suitable 
vectors. Under the above condition the system (o.2.,l) 
admits a trivial solution, which is identically equal 
to zero, we call it as null solution. 

The null solution, under certain control function u, 

is said to be 
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(i) g— sizable, if iliere exists a 6 > o such, tlial; 
for .any ||xCtQ)|j < 5, the solution of (5,2.1) satisfies 

1 I4t)t 1 < g(t) 

for all t ^ for some given non-negative continuous 
function g(t). 

Clearly the above definition coincides v/ith that of 
Lyapunov, if the above is true for any constant 
g(t) = e > 0 . 

(ii) quasi-stable if for anj?' given G > 0 , there 

exist 6 > o and some > t^ > 0 stich that ! 11 e 

for aJ.1 t > an^ 1 {^eCto)! 1 1 6 . 

!Ehe above definition can be considered as g-stable 
with g given, by 

g(t) = r "^0 1 1 '*'1 " 

= r - — ^(t-t^+r),t^-w < t < t^ 

= G > t q 

where r is some admissible upper bound and w any 
real number greater than zero. 

(iii) asympt oticall3" stable if it is g-stable 
v/ith some g as 3 rmptotic to zero. 
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(iv) e:'5)onentially assnnpt otic ally stalle if there 
exist-., a > 0, p > o such that the system is g-stable 
with g given by 


g(t) 


-aCt-t^) 

^ e 


•fc > to 


(■^) ^ stable if H 6 > o such that for any 

1 |x (to)li 5-6 the solution satisfies 


1 ix(t)l 1 <, g(t) for T t e [ to, T ] 

for some non-negative continuous function g. 

In most of the practical problems it 
is sufficient to have gg, stability, when T is large . 
In some cases to avoid mathematical complicacies we 
will only study stability for arbitrarily large I 

but not equal to infinitj’’. 

In this chapter vre use the folio v\fing 
norms, for any element x G we denote the norm by 

n p l/2 

ibii = [.fjnr] 

XI 

i.e. the Buclidean norm, for any element 0 G C , we 
define the norm by 

110 lic = 11 ^^^^^11 

. "b 
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provided the 3upreriiui!i exists, ind let the norm on the 
op;.r;itors A(t) ^md G(t,s) etc he the one generated hy 
the corresponding Euclidean norm on Domain and range 
npaceo. We represent these norms as before h 3 r I U| I 
and I j Gj I ^ respectively. 

GIOS SD LOOP ST.ailII!Y 5^3 

Let ’^'(tjt^) be the fundamental matrix of the linear 
system 

x(t) = l».(t) x(t) (5.3.1) 

with « (toyt^) = I, to = 0 , where l(t) is an n x n 
matrix v/ith intregahle elements, let 5[ i.(t) ] denote 
the largest eigen value of the matrix 

(l/2) ■li'('t) + ("^ ) ] 

Then v/e have the follov/ing lemma [ 32 ] . 

IHEU 5.3 .2 

11 , 

Let the norm for vectors in R . be the Euclidean 
norm and that on ^’(t,to) be the one generatea by 
it, then 

1 1 ^(■fc,to)l ! < exp / d[ -(s) jds (5.3.3) 

■^0 

THEOM 5.3.4 

Let the follov/ing conditions be satisfied 

(i) There exists an men matrix n(t) with oounded 
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Ki :^'.luural^le functions as elements such, that 

S[ ^(■t) ] < -a , a > 0 

(ii) Thoro exists a function Yi 1 x D -* U 
satisfying the condition (i) and (iT) of property 
(!!) mentioned in definition 3.2.9, where I,D and 
U arc open subsets of R, and R^ respectively, 
such that the function R defined by 

P(t,x) = f(t,x(t), V(t,x)) - A(t) x(t) 

satisfies the property (H) v/ith f replaced by R 

(iii) There exist a continuous bounded non-negative 
function g(t) and a & > o such that 

0 < 6 < g(t) Ij, / g(s) e“® ds 

“ 0 

Then if Ip < a, the system (5.2.1) is g-stable, for 
miy x(tQ) = Xq G R^"^. such that ■ 

lixoll 1 S 

v/ith the control u = V(t,x). 

P;jjoof ; Assume that the conditions of the 

theorem hold and let x^ be such that | jx^j j < I • 

The system (5.2.1) is clea^^ : 

i(t) = l(t) x(t) + R(t,x(t)) 
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vvhi;!*',; P(t,x(t)) io as mentioned in the hypothesis (ii), 

••vh.ich l3i ti’i’ii is c uivalent to • 

• t 

s{t) = ‘:'(t,to) Xq + / «3(t,s) I’CsjxCs)) ds 

■^0 

Thi, c-st ability of the above IE can be studied by 
rjtiidyin^j; the fiiccd point of the operator T given by 

t 

j = T;: y(t) = '.'(t,to)x + / «(t,s)l'(s,x(s))ds 

h 

over the cl.’iss of functions defined by 

Q ={ 0 : Q e 0 ^ , I le(t)| I < g(t) } 

For any x aiid y 8 Q, we have 


llTx(t) -Ty(t)l| < / 11 ®(t,s)lLllF(s,x(s))-F(s,y(s))llds 

to 

< / Lpl lx(s) - y(s)li ds 

■^0 

Hciico 

t -a(t-s) hji 

1 iTx - iy| 1 0 < I'pl b - yl 1 0 { ® - 5--I h-yl I c 

which chows, since (Ij/cc) < 1. -tiia-t I is oohtraotlon. 

iJLso ■ . ' t ^ ^ ^ : 

llTx(t)l| < 11 ?(t,tp|l ll^oll + { II *^^’"’11 ' ' 

° X 1 1f(s5x(s)| I ds 
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•at 


b + J c 

0 


-a(t-G) 


Ijil tx(s) I |ds 


< > 6 + Ij, o-at / ,„s 

0 

3\\i\}titi\tuig the c^reission for 6 from condition (iii) 

t 

i |T3c(t)| j < g(t) 

t iiUf. T.: C Q. So we sou that T satisfies all the 
contlition.j of B;'nach’3 Fixed point theorem 2. 2.1, over 
tli. 'Wit «nd liencG there exist s a 0 £ Q such that 
for the control function u = V (t, 0(t)) the system 
(5.8.1) is g-stablo for aixy initial condition x(t|^)= 
‘''At'u iIXoll < 6 . Q.g.D. 


COROlLu aRY 5.3.5 

Under the hypotheses of the theorem 5,3.4 the system 
is stable in the souse of lyapunov. 


I/e only need to show that with 
, v/e can find out a 6 ^ o such that the 


Pro of i 
g(t) = e > o 
condition (iii) v/ill be sa,tisfiod. ITovj' 

,:(t) f S(0 6“" ds = e Lp l\ e^^ds 

n 

L 


° 1 ° 
at - ' -r u F 1 „a-fc 


=e e^^-dj/a) + (lip /a)e - e[ 1-- 1 +(Va)G 
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fhe above showa that the condition (iii) will be 
satisfied if we tales any 6 > o satisfying 6^6- 

Q «E «1 } t 


gORgijLJm 5,3.6 

Under the assuii5)tions of the theorem 5.3,4 the 
system (5.2.1) is exponentially stable. 

Proofs Ue proceed along the lines of the proof 

of Corollary 5.3,5, Since o <_ (l/a) 1 there exists 

a § > o jsuch that o < :i/(a-§) < 1. let ^ be any 
such number and 

g(t) = r e"'^'^ 
for any r > o. Then 

g(t) e*^"^ - Lp / g(s)e“® ds = r l^r / ds 

0 0 


- hi [ -1]= ^ + e 

a - ^ -I a-^ 

The above v/ill be satisfied if 5 < r . 


(a-p)t, 


u-4 ^ 

Sm? »I) , 


B. .SIG BESUIiTS 5.4 

In tliis section we present few mathematical preliminari- 
es that will be required in the next section,, details 
of these materials can be found in { 3 ] 
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DEffllTITI Qg 5.4.1 

Let •;,! be a mapping of a topological space X into 
a topological space Y and let Toe a point of X. 

We say that F is lower semicontinuous (l.s.c.) at 
^o each open set G meeting rx^ there is 

a neighbourhood U(Xq) such that 

X e u(x^) .r xO G ^ 0 . 

LEPIITIflO N 5,4.2 

The mapping r is upper semicontinuous (u.s.c.) 
at x^ if for ecnh closed set G containing 
there exists a closed neighbourhood U(Xq) such that 

X G U(Xq) FxoG 

The above can be rewritten also as follows; Eor any 
given G > o, there exists a & > o such that r:xc.[ 13 Cq]0 

whenever X G [ x^^ v/here [ A denotes the closed 

r-neighb our hood of the set A. 

Yie say that the mapping r is u.s.c. 
(l.s.c.) in X if it is u.s.c. (l.s.c.) at each point 
of X G X. F is continuotis if it is both u.s.c, and 

' . DEEIHITIOI 5.4.3 

For any sot B the lower inverse of ^ is defined 
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JTB = {X : X e Z, r X A B 0 } 

Similai-ly' 1:110 upper inverse is given by 
r"^ B = { X : X e z, r X c: B } 

THEORM 5,4,4. 

necessary and sufficient condition for r to be 
u.s.c. is that the set P x is compact for each x 
and for each open set B in Y the set r‘''B is open. 

BBPIIITIOH 5.4.5 

The mapping P is said to be closed if x^ x^^, 

^n yj;j_ 6 P Xj^ implies that 

Yo e P Xo . 

T HBORM 5.4.6 

Every u.s.c. mapping is closed. 

Proof : Let 6 > o be given. Then by u.s.c. 

of P there exists a 6 > o such that for ary x G [ 

P xc fP x„ 1 . ITov/ since 'x^ x for all sufficiently 

^ 0 0 n 0 

large n, 6 { x^ which implies P x^^cp: [ P x^ ]g. 
By theorem 5.4,4 P x^ is compact, and therefore closed. 
So, since e is arbitrary and P x^ is closed, T x^v P x^ 
for all sufficiently large n. But y^^ G ^ ^n 
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e r Xq . Since — y^, -the atove implies 3 ^q 6 P x^ . 

LMU 5.4>? 

( 'Ascoli ) 

^iny set of uniformly "bounded, equi- continuous 
functions defined over a compact interva.! is compact. 

By definition, equi-continuity of a set of functions 
means that, for any given e > o if there exists a 6 > o 
such that, for all f in the set following holds 

I I f(t) - f (^') I 1 ^ S whenever [i - 6 

LMh 5,4,8 

(Pillippov[ 33 ] ) 

let the vecter function f(t,ii) be continuous, let 
the set Q(t) be closed, bounded and u.s.c, ?/ith 
respect to inclusion in t, let the vector f(t,u) 
describe a set R(t) when the vector u describes 
the set Q(t).and let y(t) be a measurable vector 
function such that y(t) G R(t) . Ihon there exists 
a measurable function u(t) such that f(t,u(t))=:y(t) 
for almost all t . 

A more general version of the above lemma can be found 
in 34 ^ • Following is a generalization of Kalcutam s 
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fixed point theorem to Banach spaces given hj Bhonenhlust 
and Karlin [ 29 ] . 

THEOREM 5.4.9 

(Bhonenhlust, Karlin) 

let S he a convex closed set of a Banach space. To 
each point x of S a non-void set T x CT S is 
assumed given. If 

(i) ^ ^ j ^n ^ ^ f n e R 2 :^ imply y e T x 

(ii) the union t) P x over all x 0 S is contained in 
some sequentially compact set T, t'.'. r, 

then there exists a point x O S such that x e P x 

OPEh LOPE SIIBILITY 5.5 

Let g(t) he a hounded positive continuous function 
defined over [ o, T ] where T is any positive real 
number less than infinity. Lot B denote the class of o 
all Lipschitz continuous function defined over [ o,T ]} 
?/ith unformly hounded Lipschitz constant M, with va,lues 
in S cr D C: R^ , where D is an open set. The value 
for M ha,s been derived in what follows. Let the set B 
he defined as 


B ={ X : X 0 6, l|x(t)| j <_ g(t), t e[ o, T ] } (5.5,1) 
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LEM . 5.5.2 

Tile set B as defined in (5.5,1) is convex and compact. 

Let x,y e B, w S [ o,l ] and let 
z(t) = X¥x(t) + (I_w) y(t) , 

[|z(t) - 2 (t)lj < wll x(t)- x(t)l| + (I-W)j |y(t)~j(¥)l |- 

< w Mjt - t|+ (1-v/) Mtt - t| = M|* t - 1 1 

I'iius z is Lipscliitz continuous ?;itli the same Lipschitz 
constant M and hence z G 0 

1 lz(t)ys=l Iwx(t) + (l-w)y(t)(l < wj|x(t)ll+ (l-w)| )y(t)| I 

< wg(t) + (l-w)g(t) = g(t) forT t e [ 0,-^1 ] 

Which shows that z G B and hence B is convex. 

Because of the hypothesis on 6, for all x G B 

I I x(t) - x(t)l 1 < Mjt - t| . Thus B is a set of 
equi=continuous functions. .iUso B is uniformly bounded, 
which follows from the fact that each x G B is uniformly 
bounded by g(t). /md hence by Ascoli^s lemma 5.4.7 
B is compact. a.E.I) . 

let U be an open subset of and Q(t) Cl. L Le a- 
bounded and closed set defined for all t G [ o, T ]. 
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Iiet Q(t) be u«s.c. in t. Let 

^2 = Sup{ l|u(t)|| , u(t) e Q(t) } 
t 

Consider tiie LE 

x(t) = f(t, x{t), u(t)) (5,5.3) 

where f: [o, !E|xDxTJ-* Let A(t) be an n z n 
matrix with integrable elements such that 

9 1 A(t) ^ <_ -a, a > 0 

By o(t, t^) we denote the fundamental matrix of the 
linear system 

x(t) = A(t) x(t) 
with 5 (t^, t^) = I, t^ = 0 . 

Define as before the function E by 

E(t, x(t), u(t)) = f(t,x(t), u(t)) - A(t) x(t) 

and assume that E satisfies the property (EP) as 
mentioned in definition 4,2.1 with f replaced by E 
and [ 0, T] d I . As before the given DB (5.5.3) is 
equivalent to 

x(t) = A(t) x(t) + E(t,x(t), u(t)) 
v/hich again is equivalent to 



103 


t 

x('fc) = '5(t,tQ)XQ + S §(t,s) ^(s, cx(s),u(, 3 ))ds 

■^0 

for some initial condition xCt^) = . 

He are interested in the existence of a measurable 
function u(t) e Q(t) and & > o such that the system 
(5.5.3) is gj stable for any x^^ G with j ] < 6. 

lor any x e B, define the mapping 1 x by 

t ■ 

1 X = {y : y(t) = o(t,t(^)xQ + / $(t,s)P(s,x,u) ds, 

^o 

u(t) e Q(t)} (5.5.4) 

Thus lx defines a set of functions corresponding to a 
given element x in B, the set being generated by 
varying u(t) over Q(t) . 

Suppose that the function g(t) satisfies the relation 

0 < & 1 [ g(^) “ ~ S g(s) ds 

a ■ o 

LEMhA 5.5.5 

Under the conditions mentioned above- 1 x C" B 
for each x G B. 

p^^oof : Bet y G 1 x. lor convenience let 


us denote l(t, x(t), u(t)) by l(t) only. 
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1 |y(t)--y(t)| I < 11 $(t,to) - tQ)ll6 


+ II / §(t,s) P(s)ds - / $(t,s)I’(s)dsl 1 

■^o 

11 / 

1. li 'K'tj'fcQ) [ l6 + Ij / § (t, s)5’(s)ds 

■^o 

-/ $(t,s)F(s)ds| 1 + 11 / §(t,s) PCs) ds 

t ' 

- / o(t,s)P(s)de|-{. + 1 / It ^(=^,8) I’(s)l|dsl 

■^0 

t 

1. 1 I [6 + / I I <5(1;, s) - 5(t,s)| I . 

■^0 ' 

|lF(s)tlds+ I / II 5(t, s)| } .I'j; “.3'(s)l 1 ds | 


Since elements of ^ (tjs) are continuous functions of t 
for each fixed --s G [o,l ] , it is uniformly continuous 
and hence there exists a constant K such that 

11 Kt,s) 5(t ,s)lj < 1 It - t 1 

for all s G [ o, T ] . Also 

11 cj(t,s) II = Sr-p = 1 

s 

s<t 


Sup 

B 

s<t 


Hence we have 
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I |y(t) - 7(T)| I < 6 + / +L r2)ds 

■^o 

+ (^1^^+ lgrg)jt - tj 

< [ 1 6 + nT(l^r^+ Igrg) + (l-ii^3_+I.2rg) ] | t - 
= Mjt - ^ 

where M = N6 + (L^r^ + ^2^2^- 
Tile above implies tliat y G 6 . Also 

"tf 

1 1 1 S 1 1 0 (■tj'to) t + •/' t 1 ^(t , s) I I . j [ ^(sjxCs) ,u(s) I I ds 

■^0 

< 6 + )+ lgUu(s)[|)ds 

■^0 

< e"“‘^6 + 6““"^ / 1 g(s)ds + 

o a 

Sub si iiu ting the value of 5, we get 

1 \ji^) 1 1 < s(t) 

Hence y G B that is r x Cl B 


imniA 5.5.6 

The mapping V defined by (5.5.4) is upper semi- 


continuous. 
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I£2l 2£* 5'roni the proof of lemma 5.5.5 it is 

clear that the set P x satisfies the conditions of 
Ascoli’s lemma, and hence is compact for each x SB. 

Also hecause of the continuity of P(t,x,u) with respect 
to u, the inverse image of apjr set in P x is open in 
Q(t), therefore the upper inverse is also open. Thus hy 
theorem 5.4.4 P is u.s.c. q.b.D . 

TH320RM 5.5.7 

Let the follov/ing conditions he satisfied. 

(i) There exists an n x n matrix A(t) with integrable 
elements such that 6[ A(t) ] 5, “cc » a > o . 

(ii) The function P defined hy 

P(t ,x(t) ,u(t) ) = f(t,x(t), u(t)) - A(t) x(t) 

satisfies the property (PP) as mentioned in definition 
4.2.1. with f replaced hy P end [ o,! ] CP. I- 

(iii) There exist Q(t) c: U C. , a closed hounded, 
u.s.c. set in t with 

Sup{ ||u(t)|| , u(t) e Q(t) } = rg , 

t 

a number 6 > o and a hounded non-negative continuous 
function g(t) such that for all t G [ ] 

/ e°‘®g(s)ds . 
o 


^ 2^2 


at 


- L 


1 


0 < 6 < [ g(t) - 


oc 
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Then for any given e with j ] < 6, there 
exists a measurahle function u*(t) e Q(t) suoh that 
the system is stable. 

Suppose that the conditions of the 
theorem hold. Define the sets B and r as in (5,5.1) 
and (5.5,4) . vTe show that B and r satisfy the theorem 

5.4.9. 

By lemma 5.5.6 T is u.s.c, and hence by theorem 5.4,6, 
it is closed i.e. r satisfies the first condition of 
theorem 5.4.9 . Since for all x6B,rxciB, as has 
been shown in lemma 5.5.5, U V x for all x is also 
contained in B wliich is compact. Thus T satisfies 
theorem 5.4,9 and hence there exists an x G B such that 

X e r X 

Since Q(t) is closed, bounded and u.s.c, and D x is its 
image, by BillippovVs lemma 5.4.8 corresponding to 
x*G r X* there exists a measurable function u''(t) 6 Q(t) 
such that X* = r (x*, u’*^) that is 

x^ = $(t,t^)x^ + / $(t,s)B(s,x*(s),u*(s))ds 

U U 1 

■to 

ilnd since x* 6 B =$► |jx*(t)| j < g(t) for all t g[ o,T 
the system 5 . 2 .1 is st able . 
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